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BURROUGHS TRUTH FUNCTION EVALUATOR 








An Analysis of a Logical Machine Using 
Parenthesis-Free Notation 


1. Introduction. LUKASIEWICz’ parenthesis-free notation’ permits a 
simple and easily mechanizable process of truth-table computation. We shall 
describe this process and prove some relevant theorems. A 10-variable relay 
machine employing this method for a two-valued logic has been constructed 
and used by the Burroughs Corporation (see frontispiece) and it is feasible 
in the present state of the computer art to construct a 25-variable machine 
which would reduce all evaluants of a 250-character formula in a two-valued 
logic with monadic and dyadic operators in about an hour. The operation of 
such a machine will be described in terms of the manipulations which it per- 
forms upon a formal propositional language. 

2. The Language. We consider any language LZ based on a quadruple 
(C, W, P, F) whose terms will be characterized in this section’. C is a finite 
non-empty set of characters. 

(Example 1: The theory will be illustrated from time to time by the use 
of a specific sample L whose characters (members of C) are the following: 
the functors N (negation), K (dyadic conjunction), A (dyadic alternation) ; 
the propositional variables, p, g, 7; the truth-constants 1 (truth), 0 (falsity).) 

Definition 1: Any finite sequence of characters (including the null se- 
quence) is a formula. 

‘A’ will designate the null formula. Except for ‘A’ and ‘x’ (to be introduced 
later), all lower case and upper case Greek letters (with and without natural 
number subscripts) will range over all characters and all formulas respec- 
tively. Juxtaposition in the syntax language will be used to denote juxtaposi- 
tion in L. We will use the following terminology pertaining to formulas. 

Definition 2: 

A (length): Z(A) is the number of character tokens in A. 

Let A= 96Y, then 

B (tail): 7;(4) = VW where i S L(A) and L(W¥) = i, or i > L(A) and 
Y = A; 

C (head): H;(A) = ® where t S L(A) and L(#) = i, or 1 > L(A) and 
=A. 


W is a function assigning to each character (element of C) an integral 
weight S 1. 

(Example 2: N has weight 0, K and A have weight —1, the variables and 
truth-constants have weight 1.) 

The following concepts are defined in terms of weight. 

Definition 3: 

A (character degree): D(é) = 1 — W(8). 

B (formula weight): W(A) is the sum of the weights of the character 
tokens of A; W(A) = 0. 

C (max. weight, tail): Wax(A) = Max(W[T7;(A)]) fori > 0. 
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D (min. weight, tail): Wain(A) = Min(W[7,(A)]) for 7 > 0. 

E (positive formula): A is a positive formula if and only if Wain(A) > 0. 

F (well-formed formula): A is a well-formed formula if and only if A is a 
positive formula and W(A) = 1. 

(Example 3: pKqN is not a positive formula; ApgNr is positive but not 
well-formed ; KpNq is well-formed.) 

Note that the non-recursive definition 3F shows that a single scan of a 
formula from right to left is sufficient to determine whether it is well-formed. 

P is a non-empty subset of C whose elements are truth-constants. ‘x’ with 
or without natural number subscripts will range over this set. P must satisfy 
the condition 


(1) W(x) = 1. 


F is a function which assigns to each character of degree > 0 a “truth 
function” ; i.e., which satisfies the conditions 


(2a) If D(é) > 0, then F(drpqy---m1)eP; 
(2b) F(x) =z. 


(Example 4: In the sample language P consists of 0 and 1; the function 
F applied to K gives F(K00) = 0, F(K01) = 0, F(K10) = 0, F(K11) = 1.) 

We can now define the ordinary concept of (propositional) variable. 

Definition 4: 6 is a variable if and only if W(é) = 1 and 4 is not in P. 

An important relation between positive formulas and well-formed 
formulas which we will need is given by: 

THEOREM I: (A) Ais a positive formula tf and only if A can be partitioned 
(i.e., divided into a sequence of disjoint segments which exhaust it) into exactly 
W(A) = 1 well-formed formulas. (B) There is at most one partition of a posi- 
tive formula A into well-formed formulas.* 

Proof: (A) The “if” part is obvious. The ‘‘only if” part is proved as 
follows. Let I be the shortest head of positive weight of a positive formula 
A. Hence, for anyi > 0, WLT;(T)] > 0; and since L is not null it is a positive 
formula. Again, since T is the shortest head of A of positive weight, WH. 
([)] S 0 (wherea = L(T) — 1); andsince TI isa positive formula, WLT:(T) ] 
= 1;hence W(I) S 1. But then W(T) = 1 and [ isa well-formed formula. 
It follows that W[7,(A)] = W(A) — 1, for 6 = L(A) — L(T); then if 7,(A) 
is not null it is a positive formula and this process may be repeated until A 
is partitioned into W(A) well-formed formulas. 

(B) The proof is by induction on L(A). If L(A) = 1, A is a single char- 
acter and ‘‘B”’ obviously holds. Suppose ‘‘B”’ holds for all positive formulas 
of length < ;consider a positive formula A with L(A) = n. If there are two 
partitions of A, one will have a first (leftmost) well-formed formula © of 
length S the length of the first well-formed formula ®¥ in the other parti- 
tion. But W(@v) — W(%) = W(Y) = Osince W(@V) = W(4) = 1. There- 
fore ¥ = A, else Wuin(@¥) = 0 implying $Y is not well-formed. If 6 = 6¥ 
= A, then “B” has been established for A. Otherwise A = $F = SVT where 
T is a positive formula and L(T’) < L(A). In this case the inductive hypothe- 
sis shows that “‘B” holds for T and hence for © = A. 

Theorem I makes possible the following definition. 
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Definition 5 (partition function): If A = A,z---A;---A; and each A; is 
well-formed, then P;(A) = Aj. 

(Example 5: For A = pKpqNr, P:(A) = p, P2(A) = Kpgq, P3(A) = Nr.) 

With the aid of Theorem I and Definition 3A it can be shown that our 
formulation of the language L is easily reducible to a more conventional 
formulation. For example, A is a well-formed formula if and only if it is of 
the form 6Ap,)---A:, where each Ag is well-formed. Note that D(é) is the 
number of well-formed formulas following 6 in the above decomposition ; if 
D(é) > 0, 6 is an operator of the propositional calculus and D(6) is its degree 
in the ordinary sense, while if D(é) = 0, 6 is well-formed by itself and hence 
is a truth-constant or propositional variable. 

An interesting property peculiar to the parenthesis-free notation is 
given by: 

THEOREM II: Every formula of a language L is a segment of some well- 
formed formula of L if and only if L contains at least one character of negative 
weight. 

Proof: The proof of the “‘if’’ part is as follows. Let w be a character of 
negative weight. For A an arbitrary formula of the language L, VA¢® is well- 
formed where ® consists of 1— Wmin(A) occurrences of x and ¥ consists of a 
formula of weight —1 (e.g., w followed by —1— W(w) occurrences of x) re- 
peated W(A®) — 1 times. The proof of the “only if” part is as follows. 
W(x) = 2, and, since rm is a segment of some well-formed formula A, L 
must contain a character of negative weight in order that W{A) = 1. 

The need for a machine to do truth-table computation occurs only for 
languages containing at least one symbol of negative weight and a multi- 
plicity of symbols of weight 1. 

3. The Machine. With the aid of the two following definitions we can 
describe, for any given language L, the construction of a machine to evaluate 
formulas of L. 

Definition 6 (the set of specification functions) : G is a specification func- 
tion if and only if (1) if 6 is a variable, G(6)eP and (2) if 6 is not a variable, 
G(6) = 6. 

(Hereafter ‘S’ will represent an arbitrary specification function.) 

Definition 7 (evaluation function) : Es(A) = A; if 6A is a positive formula, 
Es(6A) = F(S(6)Hpe[Es(A) ])TaLEs(A) ] where a = L[Es(A)] — D(d). 

(Example 6: For the formula of Example 5, A = pKpgNr, the recursive 
evaluation is as follows: let the specification, S, of the variables be p = 1, 
q=0, r = 0, then Es[7,(A)] = FLS(r)] = 0, Esl72(A)] = F(NO) = 1, 
-++, Es[7;(A) ] = F(K10)1 = 01, Es(A) = F(p)01 = 101.) 

The following lemma shows that L(Hp )[Es(A)]) = D(é) and hence 
that Es is defined for all positive formulas. 

Lemma: If A is a positive formula, then E.s(A) is of the form m,m,_1-° +71 
where » = W(A). 

Proof: That all characters are x’s follows directly from the conditions 
(2a and 2b) which characterize F. That » = W(A) is established by an in- 
duction on L(A). 

Note that if A is well-formed, L[Es(A)] = 1; it is easy to show that 
when A is well-formed, Es(A) is its truth-value relative to S in the ordinary 
sense. 
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A number of well-formed formulas may be evaluated concurrently by 
juxtaposing them to form (by Theorem I) a positive formula A. Let A = 
5z.a)* + *5;-- +51, let m be the number of truth-constants in the language L 
(i.e., L is an m-valued logic), and let v be the number of distinct variables of 
A. The machine makes m” scans of A, one for each distinct specification of 
the variables, producing each time Es(A);a single complete scan for a given 
S with the accompanying computation determines the “Sth” major cycle. 
Each major cycle is divided into L(A) minor cycles, the ith minor cycle en- 
compassing the processing of the ith character 6,. 

The machine consists of two basic components, a Memory and an 
Evaluator. The Memory (e.g., a magnetic drum, an acoustic delay line) 
stores A and during each major cycle sends it characters 6,,4)---6; to the 
Evaluator in order of ascending subscripts. During the Sth major cycle the 
Evaluator realizes the recursive function E.s(A) by producing successively the 
E;(.7,(A) ]. It does this by means of three parts: a Specifier, a Function 
Switch, and a Register. 

The Specifier (e.g., an electronic counter with switching gates) effects the 
sequence of S functions required by A; during the ith minor cycle it receives 
5; and produces S(é;), and at the end of the Sth major cycle it introduces a 
new S. During the 7th minor cycle of the Sth major cycle the Function Switch 
(e.g., an array of electronic switching gates) receives S(6,;) from the Specifier, 
Ape; (EsT:-1(A) ]) from the Register, and produces the new character 
F(US(6;)Hpe;)(Esl7:-1(A) ]) ]. This new character is sent to the Register 
(e.g., an electronic shifting register) which by shifting an amount W(é,) (a 
positive value indicates a right shift, a negative value a left shift) substitutes 
it for Hp@;)(EslTi-1(A) ]}). 

4. Some Theorems. The following theorem justifies the use of a positive 
formula A to compute the truth-values of its well-formed components P; (4). 
THEOREM III: Jf A is a positive formula then Es(A) = Es{Pw)(A)}:- 

Es[P,(A) ]. 

(Example 7: In the formula of Examples 5 and 6, Es(A) = Es(p)Es 
(Kpq)Es(Nr) which, for the S of Example 6, evaluates to 101.) 

Proof: The theorem follows directly from the fact that, if @ and W are 
positive formulas or null, £s(@V) = Es(#)Es(¥). This may be proved by an 
induction on L(#). It is obviously true for 6 = A. We assume it to be true 
for some [ and show that it holds for yl’. By Definition 7 Es(yTWV) = F(S(y) 
Ap lEs(T¥)]) (Tal Es(T'¥) ]), where a = LLEs(f'¥)] — D(z), and by the 
inductive hypothesis Es(yP'W) = F(S(y)HpwylEs(M)Es(¥) )) (MLEs(T)Es 
(¥)]) where 6 = LLEs(T)Es(¥) ] — D(y). Since 71 is a positive formula, 
L{Es(I)] = D(y) and we have Es(yP¥) = F(S(y)HoqLEs(l) DHLEs 
(T) JEs(W), where c = L[Es(T)] — D(y); hence Es(yF'¥) = Es(yV)Es(¥). 

For a given positive formula A the Register must be able to store a 


(4) 

(positive) formula of length Max [L(Es[T;(A) ])]. The following two 
theorems show how this quantity —— upon the. structure of A. 

THEOREM IV: Jf A is a positive formula then Max (L(Es(T;(A)])] = 
Wwax(A). 

L(A) =6 
(Example 8: In Example 6 it can be seen that Max [L(Es[7;(4)])] 
i=1 
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= L[Es(A)] = 3 = Wwax(A). It may be checked that the evaluants of 
T;(A) and 7,(A) which were not given will not change this maximum.) 
Proof: By the Lemma and Definition 3C. 


WA) 
THEOREM V: If A is a@ positive formula, then Wmax(A) = Max (Wwsx 
j=i 


[P;(4)] +7 — 1). 

Proof: Consider all the tails of A and note that each complete P;(A) in a 
tail contributes a weight of 1. 

This theorem shows that, given a set of well-formed formulas to be 
evaluated as a positive formula or to be used as arguments for a commuta- 
tive operator, a formula of smallest Waa. may be formed by placing these 
well-formed formulas in order of ascending Wmax from left to right. 

In designing a specific machine to carry out our process some decisions 
must be made concerning the relative sizes of the Memory and the Register. 
Our concluding theorem gives information relevant to this decision. 

THEOREM VI: Let My.p be the set of well-formed formulas V such that (1) 
the maximum degree of any character of V is D > 1 and (2) Wmax(¥) = W. 
Then (A) for any V «¢ Mw.p there is a ® € Mw.p such that L(®) = L(¥) and 


(B) L(®) =W+ - et where {z} is the smallest integer = s. 





D-1}’ 

Proof: For part (A) consider for any ¥ « My,p the formula ¥; = dxp--- 
m2¥, where D(5) = D. Clearly V; satisfies (1). By Theorem V, Wmuax(¥1) = 
Max(D,W) which equals W since, for any positive formula A, Waax(A) = 
the maximum degree of any character of A. Hence ¥; ¢ Mw, p and L(¥;) > 
L(¥). To establish (B) consider any VW ¢« My,p. By (2) Y must contain at 
least W characters of weight 1. Since W(¥) = 1 the sum of the weights of 
its characters of negative weight is at most 1 — W; by (1) no character of 
W has a weight less than 1 — D; since D > 1, YW must contain at least 








W-1 he 
| Dn1 | characters of negative weight. Hence L(¥W) = W + D-1 } 


For a language L containing characters of all degrees between 2 and D 
inclusive, if Mw,p is non-empty, it contains a well-formed formula of the 


A : : W-1 . . 
minimum length which consists of Du-1 characters of negative weight 


followed by W characters of weight 1. 
ARTHUR W. BurRKS 
Don W. WARREN 


Jesse B. WRIGHT 
University of Michigan 
Ann Arbor 


The writing of this paper and the research which it reports were done under the sponsor- 
ship of the Burroughs Corporation. 

1 JAN LUKASIEWICZ, Aristotle's Syllogistic from the standpoint of Modern Formal Logic. 
Oxford, 1951, p. 78. 

2 In this section we employ the results of STANISLAW JASKOWSKI, KARL MENGER, KARL 
Scur6TeR, and D. C. GeERNETH. See Paut ROSENBLOOM, The Elements of Mathematical 
Logic. New York, 1950, ch. 4, sec. 1, for a discussion and further references. 

* Theorem I implies that under the operation juxtaposition the set of positive formulas 
together with the two-sided identity A is a semigroup which is generated by the set of well- 
formed formulas. 

_ 4For a description of computer components see ENGINEERING RESEARCH ASSOCIATES, 
High Speed Computing Devices. New York, 1950. 
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Polynomial-Like Approximation 


Frequently, in machine computations, one may substitute for a function 
of relatively complicated analytic structure a function of simpler form with- 
out appreciable error. Moreover, if the function is given tabularly, a con- 
siderable reduction of tabular input data (load on the memory) may be 
achieved. In what follows, we deal with real-valued functions of continuous 
real variables x, y,---. Computationally, it is usually better to work with 
variables assuming finitely many discrete values. In the formulas which fol- 
low, however, the discrete analogues are apparent, integration to be replaced 
by summation, etc. We have in mind multivariate approximations by poly- 
nomial-like forms, in particular, by functions M of the form 


(1) M(x,y) = E d1(e)¥s(0), 


in the bivariate case. Some of the following results are generalizable to func- 
tions of more than two variables. 

Of special interest is the so-called slide rule form, f(x) + g(y). We shall 
give formulas for approximating to a given real-valued continuous function, 
defined over the unit square 0 < x, y < 1 to minimize, respectively, each 
of the following measures of goodness of fit: 


(a) f f [s(x, y) — f(x) — g(y) Pdxdy, and 
* max, |2(e, 9) — £2) 80) 


A simple variational technique yields for (a): 
1 
f(z) = f s(2, vay + co, and 


1 
g(y) = f 2(x, y)\dx + ¢:, where 


Co and c; are any real numbers such that 


1 1 
ota=- f f 2(x, y)dxdy. 


Recursive formulas for (b) are based heuristically on the fact that given 
two numbers a and 3, one can approximate to them by a single number c; 
taking c to be the arithmetic mean of a and 6 minimizes the maximum abso- 
lute error. The recursive formulas for (b) are as follows: 


fnti(x) = § { max [2(x, y) — ga(y)] + min [2(x, y) — ga(y) J}, 
gn(y) = § { max [2(x, y) — fa(x)] + min [2(x, y) — fa(x) J}, 


and fo arbitrary but continuous, ” = 0, 1,---. Dir1pertTo & Straus! have 
shown quite generally that the above process converges to a pair of continu- 
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ous functions (f, g) which has the required property. Included in their dis- 
cussion is a precise method, discovered independently by the writer, for 
estimating the value of the minimum. In this connection, the reader is 
referred to their paper. 

We now return to the more general form (1). Let » > 0 be a prescribed 
integer. We assume that z (continuous over the unit square) is not of the 
form (1) for this particular ”. We shall give two sets of formulas for approxi- 
mating to z by continuous functions of the form (1), which satisfy respec- 
tively: 

(c) The error vanishes on some rectangular grid of lines x = x1, ¥2,°--, Xa, 
Y = V1, V2,°**, Wn (x; F x; and y; ¥ y; if 4 ¥ 7) 


(d) f f [s«. y) - Fs o,(e ws) | dedy is minimized. 


For (c), functions ¢;, ¥; can be computed successively as follows. Choose a 
point (x1, yi) such that 2(x,, y:) * 0 and two constants ¢, ¢;’ such that 
C101’ = 1/2(x1, y1) and set 


oi(x) = c12(x, yi) and vi(y) = cy'2(x1, y). 


We compute ¢2 and yz in precisely the same manner, except that z is now 
replaced by a new function 22 given by 


Zo(x, y) = 2(x, y) — oi(x)~r(y). 


Thus, we choose a point (x2, y2) such that z2(x2, ye) # 0, and two constants 
C2, Ce’ such that cece’ = 1/22(x2, y2) and set 


g2(x) = Coze(x, y2) and 
v2(y) = C2's2(x2, y). 
We continue in this manner, with 
Ser1(x, ¥) = 2u(x, y) — de(x)ve(y) 


until we have the required number of functions. A simple inductive argu- 
ment shows that condition (c) is satisfied. 

That the solution of problem (d), i.e. finding a least squares fit to z of 
the form (1), bears a direct relationship to the theory of integral operators 
in Hilbert space was pointed out by G. W. Brown a few years ago. His 
observations are essentially these: 

Let Ax, A2,- ++, An be the largest n eigenvalues of the symmetric kernel 


K(x, y) = f 2(x, t)2(y, t)dt 


and ¢:, ¢2,°--, $n, be corresponding eigenfunctions (normalized so that 


1 
f 7 (x)dx = 1) and set 
0 


vio) = fst, ndos(a)de, jf = 1,2,0--49, 
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then the @’s and y’s defined in this way form a solution to (d), the ¢’s being 
mutually orthogonal and the y’s being mutually orthogonal. The computa- 
tions can be carried on sequentially. Let A, be the largest eigenvalue of K, ¢; 
a corresponding normalized eigenfunction and y; defined as above. Having 
found ¢; and ¥ we form the symmetric kernel, K; associated with the 
residual function 2(x, y) — $1(x)¥i(y). Then, the largest eigenvalue of K; is 
the next largest eigenvalue of K and we take ¢2 to be the corresponding 
normalized eigenfunction with y. defined as above. We continue in this 
manner until we have the required number of functions. Incidentally, for 
n = 1, a variational technique yields the necessary conditions for an 
extremum: 


oe) = fst, dworas/ f" vores, 


vo) = fae, no(erdr/ f" eee. 


These can be used to generate an iterative computation for ¢ and ¥. How- 
ever, questions of convergence, proper initiating functions to achieve the 
largest eigenvalue, etc. seem to be difficult. One final remark—the minimum 


1 1 n 
value of f f [z(x, y) — ¥ $;(x)¥;(y) Pdxdy is precisely the sum of the 
0 7=1 


remaining eigenvalues of K, i.e. the sum of the eigenvalues minus the sum of 

the largest ” eigenvalues, the former sum being equal to { /2*dxdy. 
OLIVER Gross 

The Ranp Corporation 

Santa Monica, California 


1S. P. DitiBErto & E. G. Straus, “On the approximation of a function of several vari- 
ables by the sum of functions of fewer variables,”’ Pacific Jn. of Math., v. 1, p. 195-210, 1951. 


A Logarithm Algorithm 


The method of calculating logarithms given in this paper is quite unlike 
anything previously known to the author and seems worth recording be- 
cause of its mathematical beauty and its adaptability to high-speed com- 
puting machines. Although there are well known methods! which involve 
continued fractions, these methods invariably utilize the analytic properties 
of the logarithm function and not the arithmetic properties of the individual 
logarithm. The first version of this algorithm is based directly upon such 
arithmetic continued fractions. In a subsequent skeletonized modification, 
however, continued fractions no longer appear explicitly. 

Let ao > a: > 1 be given. To find log,,a; we determine the two se- 
quences 

ado, a3, ee 
M1, N2,°**, 


where the n’s are positive integers, by the relations 
a** < a1 < a," 


Qin, = Ay1/a,™. 


(1) 
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eing We define the complete quotient x; > 1 by 
uta- niti/z1 
s $1 Si ao = a ° 
ving | ince ee 
the . , . 
¢ is we write 
ding Xx, = M2 + 1/x2 
= and so on. Thus we have the continued fraction 
, for 
an wee oF ee 
Ba nt mat mateo 
For example let a) = 10, a; = 2. Then 
2? < 10 < 2+. 
Therefore nm; = 3 and az = 10/2* = 1.25. Further: 
low- 1.25% < 2 < 1.254. 
the 
num Thus m2 = 3 and a3; = 2/1.25* = 1.024. Continuing we obtain 
the 1 n; a; 
0 — 10 
2 3 1.25 
: 3 9 1.024 
4 2 1.009741958 
, 5 2 1.004336279 
vari- 
1951. 6 4 1.001041545 
It follows that 
= tee we 
me 3458494549444 «>. 
slike Breaking off the continued fraction at successive terms we find: 
be- 1 rational approximation to logio2 
om- 1 1/3 = .33333333 
olve 2 3/10 = .30000000 
rties | 3 28/93 = .30107527 
dual 4 59/196 = .30102041 
such 5 146/485 = .30103093 
on, 6 643/2136 = .30102996. 
. se- Since logio2 = .30102999566 we see that each successive approximation 


one. This rate of convergence, one decimal place per cycle, is typical. If one 
chooses at random an dy» and an a, one will “almost always’’ find such a 
rate of convergence. A method of estimating this universal rate is given 
below. 

Both phases of the recurrence rules given in (1) may be carried out simul- 
taneously. One divides a;_, by a;; then the quotient by a;; then that quotient 


| gives us approximately one more correct decimal place than the previous 
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by a,, etc., until the resulting quotient is less than a;. Then this last quotient 
is @;,, and the number of divisions is m;. In fact, this part of the process in- 
volves only division. 

The second part of the process, namely, that of obtaining the successive 
fractions, may similarly be reduced to only addition, but at first we note 
that from the theory of continued fractions if P;1/Q;-1 and P;/Q; are the 
rational approximations which include the terms up to 1/n,_; and 1/n; 
respectively ; then 


(2) Piss Pia t+ mis 


Qi41 " Oia + ins 


For example, above we had 


P2/Q2 = 3/10, P3/Q3 = 28/93, m = z3 





and therefore 
Ps/Qs = 59/196. 


However, for an automatic computing machine we recommend the fol- 
lowing variation which abstracts the contents of the last two paragraphs. 
We maintain six registers A, B, C, D, E, and F. At each inning we do one of 
two things: 


Operation I (if A > B): 

We put A/B in A, C+ Ein C,and D+ Fin D. 
Operation II (if A < B): 

We interchange A and B, Cand E, D and F. We start with ap in A, a, in 
B, 1 in Cand F, and 0 in D and E. The latest approximation to log,,@; is 
always E/F. 

In the example above we would obtain: 


A B Cc D E F 
10 2 1 0 0 1 
5 2 1 1 0 1 
Op. I 2.5 2 1 2 0 1 
. 1.25 2 1 3 0 1 
Op. II — 
2 1.25 0 1 1 3 
Op. I 1.6 1.25 1 4 1 3 
; 1.28 1.25 2 7 1 3 
1.024 1.25 3 10 1 3 
Op. II —> 
1.25 1.024 1 3 3 10 


It is readily seen that in this variation we need not assume dp > 4, as 
was done in (1), but merely that ap > 1 and a; > 1. Further, if one or both 
of these numbers are less than 1, then since 


logax-? = — logax, logijex = — logex, logijax = logax, 


we may proceed as before after taking the reciprocals of those numbers less 
than 1. We then multiply the fraction obtained by (—1)™ where m is the 
number of reciprocals taken. 
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If it happens that the logarithm is a rational number, for instance loge4 
= 2/3, then at some point B becomes 1, the exact log is obtained and no 
further changes in E or F occurs. For example: 


NNNHeE NK PD OOM 
Ree Dh Pb by 
Owen ORE 
NPR WHR OY 
NNN KE OO fy 
AE 


If only four registers are available then one may (with some care) 
economize by keeping C and D side by side in a register C and similarly E 
and F in a register E. 

Operation I now reads: _ 

We put A/B in A and C + Ein C. 
Operation II now reads: s 

We interchange A and B, C and E£. 
The example now reads: 


A B Cc E 
10 2 -00010000 -00000001 
5 2 -00010001 -00000001 
2.5 2 -00010002 -00000001 
1.0010415 1.0043363 06432136 .01460485 
Op. Il 
1.0043363 1.0010415 -01460485 -06432136 


We now split .06432136 into two parts and divide: logio 2 ~ 0643/2136 = 
.30102996. Of course the “D” and ‘‘F’’ parts of the numbers must not be 
allowed to overlap the first halves. 

The simplicity of the rules given in either pair of Operations I, II isa 
recommendation for use of this logarithm in automatic computing machines. 
So also is the fact mentioned above that for ‘“‘almost all’’ ap and a, the rate of 
convergence is essentially independent of these numbers and is, on the aver- 
age, about one decimal place per complete cycle (that is, for each Operation 
IT). We justify this claim as follows. KHINTCHINE? has shown that for ‘‘al- 
most all’ real numbers, 0 < x < 1, the regular continued fraction has the 
property that the geometric mean of the m’s is an absolute constant. This 
number, known as Khintchine’s constant, has been computed by the author 
to ten figures*. Its numerical value is approximately 


K = 2.685452001. 


Now we saw in (2) that the ’s regulate the rate of growth of the Q’s. To 
estimate this growth we assume that the Q’s are governed (in the mean) by 
the linear difference equation : 


Qiss = Qi-1 + KQ. 
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From the theory of such equations, we find that for large 7: 
Q:; = a{3(K + [4+ K*})} 


Now it is easily seen that the error in any approximation P;,/Q; is of the 
order of (Q;)~* and therefore, per average cycle, the error should decrease by 
a factor of 4{K + (4+ K*)}~ = 1/9.1 or approximately one decimal 
place. 
DANIEL SHANKS 

Naval Ordnance Laboratory 
White Oak, Maryland 

1See for example, D. TEICHROEW, ‘Use of continued fractions in high speed com- 
puting,” MTAC, v. 6, 1952, p. 127. 

2 A. KHINTCHINE, “Zur metrische Kettenbruchtheorie,” Compositio Math., v. 3, 1936. 
p. 276-285. 


‘ DANIEL SHANKS, “‘Note on an absolute constant of Khintchine,” MTAC, v. 4, 1950, 
p. 28. 


The Product Form for the Inverse in the 
Simplex Method 


Summary: When a matrix is represented as a product of ‘‘elementary”’ 
matrices, the matrix, its transpose, its inverse and inverse transpose are 
readily available for vector multiplication. By an “elementary matrix’’ is 
meant one formed from the identity matrix by replacing one column; thus 
an elementary matrix can be compactly recorded by the subscript of the 
altered column and the values of the elements in it. In the revised simplex 
method,! both the inverse and inverse transpose of a ‘‘basic’’ matrix are 
needed; more significant, however, is the fact that each iteration replaces 
one of the columns of the basis. In the product form of representation, this 
change can be conveniently effected by multiplying the previous matrix by 
an elementary matrix; thus, only one additional column of information need 
be recorded with each iteration. This approach places relatively greater 
emphasis on “‘reading’’ operations than “‘writing’’ and thereby reduces com- 
putation time. Using the I.B.M. Card Programmed Calculator, a novel 
feature results: when the inverse matrix is needed at one stage and its trans- 
pose at another, this is achieved simply by turning over the deck of cards 
representing the inverse. 


Introduction: 


The simplex method is an algorithm for determining values for a set of 
non-negative variables which minimizes a linear form subject to m linear 
restraints.!-2¢5 It may be characterized briefly as a finite iterative procedure. 
Each iteration produces a new special solution to the restraint equations 
involving a subset of m of the variables, only one element of the subset 
changing on successive iterations ; the remaining » — m variables are equated 
to zero. The vectors of coefficients corresponding to the subset of m variables 
are linearly independent and constitute a basis in m-dimension real vector 
space. In the original simplex method”* (as coded for the SEAC*‘ or as found 
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in CHARNES et al.*, it is required that all the coefficient vectors be represented 
in terms of the latest basis; since the changes of basis are step-wise, a simple 
recursion relation suffices to alter the representations on each iteration. 

The revised simplex method!'-' differs from the original method in that it 
uses the same recursion relations to transform only the inverse of the basis 
for each iteration. It has been introduced to reduce the quantity of writing 
at each step (which it does in general), and is particularly effective for linear 
programming models where the original matrix of coefficients is largely 
composed of zeros, as for example, in the transportation model* or dynamic 
economic and production models.” If the original method is used, these zeros 
would be replaced by non-zeros in the successive iterations and this greatly 
increases the computational effort. On the other hand, the revised method 
leaves those zeros intact. 

One important feature of the simplex method is concerned with the cri- 
teria by which one of the vectors in the basis is replaced by a vector not in 
the basis to form the basis of the next iteration. When the constant terms of 
the restraint equations are not general, the choice of the vector to drop from the 
basis may be ambiguous and an arbitrary selection (as pointed out in unpub- 
lished examples by ALAN HOFFMAN and PxiLip WOLFE) may lead to non- 
convergence. Several devices exist, however, for perturbing the constant 
terms so as to avoid this difficulty. The earliest proposal along these lines** 
consisted in modifying the vector of constant terms by a specially weighted 
combination of the unit vectors. This approach may be used conveniently 
both for the revised and original simplex methods'-*. With the original sim- 
plex method, there is another natural way to form the pertubation which 
consists in adding a weighted linear combination of the column vectors to 
the vector of constant terms. This was suggested first by ORDEN and de- 
veloped independently by CHARNEs*. 

Although considerable attention has been paid to the above difficulty 
(called degeneracy), it usually does not lead to non-convergence. The type of 
problems in which it can cause non-convergence appear to be exceedingly 
rare. To date, there have been only two examples and these were artificially 
constructed for this purpose. Accordingly, the SEAC code and the RAND 
code use an arbitrary selection criteria in case of ambiguity. In these codes, 
a deliberate decision was made to use a simple code in lieu of a more complex 
one needed to cover a possible case that may mever arise in practice. 

The present method of using the product form for the representation of 
the inverse of the matrix, also makes use of this simplification. Again, pro- 
vision could be made for covering the rare non-convergent case, but again, 
it does not appear to be worth-while. 

We shall now describe a process by which a square non-singular matrix 
may be expressed as a product of elementary matrices of the form (2) below. 
This is illustratively seen for the simplex process which involves a step-wise 
change of basis matrix, that is to say, two successive matrices differ by only 
one column. Using a notation consistent with'® let BO-? = (Po, Pj,,---, 
P;,,) denote the (J — 1)* basis. If, in the next basis, P, is to replace P;,, then 
it is easy to show that 


(1) [BO}! = ELBOo}, 
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where E; and E;~' are elementary matrices related by 


(2) Ei = [Uo,: 7 a Vous Nl, Ur+1)° shes Un | 
as [Uo,- "4 U,-1, Yi, Ursiy* al Unt", 


where U; is a unit vector with unity in the (i + 1)** component, is a 
vector whose components 7;; are related to components y;; of Y; by 


nit = — Vi/Vrty 1¥#r 


3 
( ii = 1/yr2, 


where it is necessary that y,; ~ 0 and Y; is defined by 
(4) Y; = (B@»}'P,. 

Successive applications of (1) for / = k,k — 1, k — 2,---, 1 yield 
(5) (B®)! = ExEys---E(BO}, 


where B is the initial basis. It is usually easy to arrange that the initial 
basis B be the identity matrix so that B® may be dropped from (5). 
Consider the problem of computing a row vector 8, defined by 


(6) Bo = aB = ak, Ey-1--- Ei, 


where a is a given row vector (actually a unit vector’). Such a vector is 
required by the revised simplex method as the first step in determining the 
vector P, to introduce into the basis. It is clear that 8) can be obtained by 
successive transformations on row vectors, i.e., forming (a)Ex, (aEx)Ex-1, 
+++, etc. However, when a row vector A = (do, @1,---, Gm) is transformed 
into a row vector B = (bo, bi,---, bm) by multiplying A on the right by an 
elementary matrix E; one obtains simply 


b; = a, 1 # T1 
b, - Zz nitd:, 
0 
where, because r may be different for different /, we have set r = r;. 
Consider next the problem of computing Y by relation (4). 
(8) Y = B-'P, = E,E,-1°--E;P,. 


It is clear that Y can be determined by successive transformations on column 
vectors, i.e., forming E,(P,),---, etc. However, when a column vector C = 
{Co, €1,++ +, Cm} is transformed into a column vector D = {do, d,---, dm} by 
multiplying C on the left by a matrix of the special form E, one obtains 
simply 


dy = 65 + nitlry, ir 


(9) 
d,, = Nr Jc, r 
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From (7) and (9) it is clear that the only essential information contained in 
E; is the set of values ;; and the index r;. Note further that in (8), the suc- 
cessive E; are used with increasing I and it follows from (9) that it is necessary 
to know r; before using the n,. On the other hand, in (6), the Z; are used in 
decreasing sequence of / but from (7) it is not necessary to know r; until 
after the ni; have been used. The perfect complementarity of the preceding 


two sentences, together with the fact that >> 7,0; can obviously be com- 


0 
puted starting with 7 = m as well as with i = 0, makes it clear that (6) may 
be computed using the information in the reverse order of that used in (8). 
Let L; denote the ordered set of ‘words’ of information 


(10) li = {r13 nor, M1y***y Nt}. 


Then each change of a column of B will produce a new L:,; which may be 
stored in consecutive order to the previously computed Ly, L2,---, Li. 

On the CPC, by punching two sets of instructions on each card—one 
being, in form, the reflection, in the vertical center line, of the other (with 
appropriate adjustments for difference in algorithms (7) and (9))—the 
transpose use of the inverse may be accomplished by simply turning the cards 
over using the vertical center line of the card for the axis. 

GeEorGE B. DANTZIG 
Wma. OrcHArD-Hays 
Ranp Corp. 
Santa Monica, Calif. 


1G. B. Dantzic, ALEX OrDEN & Puitip WoLFE, The Generalized Simplex Method. 
Ranp P-392-1 dated August 4, 1953. 

2 T. C. Koopmans, Ed., Activity Analysis of Production and Allocation. New York, 1951. 

(a) GeorGEe B. DantziG, Maximization of a Linear Function of Variables Subject to 

Linear Inequalities. P. 339-347. 

(b) The Programming of Interdependent Activities: Mathematical Model. P. 19-32. 

(c) Application of the Simplex Method to a Transportation Problem. P. 359-373. 

(d) T. C. Koopmans & S. REITER, A Model of Transportation. P. 222-259. 

3A. CHARNES, W. W. Cooper & A. HENDERSON, An Introduction to Linear Program- 
ming. New York, 1953. 

4A. HorrMaAn, M. Mannos, D. Soxotowsky & N. WIEGMANN, “Computational 
experience in solving linear programs,” Soc. Industrial and Applied Math., Jn., v. 1, 1953, 
p. 17-33. 
- wo B. DantziG, Computational Algorithm of the Simplex Method. Ranp P-394, April 

, 1953. 


On Modified Divided Differences II 


[Continued from MTAC, v. 8, p. 1-11] 


Errors of Type (c). A question that presents itself is the extent to which 
errors of Type c will mask an isolated error. It will be desirable to approach 
the problem from a statistical standpoint, and to introduce the simplifying 
assumptions that the errors of Type (c) behave like round-off errors, subject 
to the following restrictions: 
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(a) The errors e¢; in the entries 1; all have the same, uniform distribution 
between —4a + c¢ and 4a +c, where a and ¢ are fixed constants. Thus if 
entries are rounded to a fixed number of decimal places, the assumption is 
that the rounding will range uniformly between +3 units of the last place; 
that is, c is zero. Or, if the entries are ‘‘chopped’’—that is, all digits beyond a 
fixed decimal place are dropped, without rounding, then the assumption is 
that the error in the last place ranges between 0 and unity; that is, c = }. 
It will be shown that the distribution function is independent of c, provided 
c is the same for all 1,. 

(b) The errors in the tabular entries are independent of one another. 
The conditions (a) and (b) constitute a useful idealized model. 


For the case of equally-spaced arguments the distribution function of the 
round-off error in differences of the first, second, and third order has been 
given explicitly by Lowan & LADERMAN!, and the method can be used for 
obtaining the distribution function for differences of all orders. A somewhat 
more elaborate study has been published by A. VAN WIJNGAARDEN?’. We shall 
here follow the method of Lowan and Laderman, based on Fourier trans- 
forms. Consider the sum 


(2.16) Wn = Aco + Aiei +°--+ An€n, 


where the coefficients A; are constants, and all the values of e, are subject to 
the restrictions (a) and (b). Let f(w, x)dx denote the probability element of 


t 
the distribution ; that is, f f(w,x)dx = F(w,t) is the distribution function 


of w. For the case when c = 0 in condition (a), we have 


Rs 1 1 : 1 
(2.17) fle,x) = a if — 7 <x< 92: flex, x) = 0, if |x| > 32 


a 


(2.18) f(Aes, x) = 7A 


, 1 1 
, if — qalA| <x< RtlAl, 
1 
f(Aex, x) = 0, |x| > aA, for constant A. 


The characteristic function g(w, ¢) associated with a distribution function 
F(w, t) is defined by 


(2.19) g(w,t) = fe fw, »)dx, 
and by the Fourier inversion theorem 

1 Pris 
(2.20) f(w, x) = xf e~* o(w, t)dt. 


It is known that the characteristic function associated with the distribu- 
tion of the sum of m independent variables is the product of the characteristic 
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functions associated with the distributions of the m individual variables. 
This gives, for w = Ae, 


$o1 4! tiie as sin (4aAt) 
(2.21) g(w, t) = Fa baa a 
Hence for w, defined in i 16) 
(2.22) yin. 0 < fae 


k=0 3aA xt 
Using (2.20), the frequency function for w, is given by 


i 1 _ sin (3aA;£) 
(2.23) f (Wa, x) = ta t ee 


2x 
1 (cos (tx) * sin (4aA,t) 
- if I 2 


T {*+ 1 k=Q 4aA k 


dt 





dt. 


The probability that w, takes on a value between 6 and c, c > 3, is then 
given by 


f f (Wa, x)dx. 


It is to be noted that the integrand in (2.23) is unchanged when A; is 
replaced by —Aj;; it will therefore be convenient to write 


Bt= aA, , 
eae G(x, t) = cos tx- ee 
bmp 2D 
Let 
G® (x, t) = d'G(x, t)/dt*; G(x, t) = G(x, é). 
Then 


G® (x,t) = t!*V(x,t), OS Ron+1, 


where V(x, ¢) is bounded for all ¢«. Moreover G® (x, ©) is bounded. Integrate 
(2.23) by parts. This gives 


G t 1 @ GW) (x, t 
(2.25) f(we, x) = +| ai +if oe at 


The term within the bracket vanishes at both limits and only the integral 
remains on the right-hand side of (2.25). Repeat the process of partial in- 
tegration until we arrive at 


(2.26) Sine tie 1° s 1 G(x, GO (x, t) 


n' t 


By expanding G(x, #) into a sum of sines (or cosines) and performing the 
ee with respect to #, it can be verified that f(w,, x) = 0 if |x| 


£= 9 =F 1Aul, as it should. In the case where the numbers A, are the binomial 
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coefficients (— 1)* + 


ence of u,, due to the individual errors e,. In the case of modified divided 
differences, the numbers A; represent M,,:;, in even central differences or 
corresponding coefficients in differences of odd order. It should be observed 
that in all divided differences, whether modified or not, 


), W, represents the error in the wth ordinary differ- 


(2.27) Mi, = 0. 


For the coefficients M;,, are independent of the function ~. In the special 
case when 4, = 1, the differences must vanish; hence (2.27). As a conse- 
quence of (2.27) it is clear that if e, = c + 6, where 6; is uniformly distrib- 
uted with mean zero, the constant c drops out in the sum represented by w,. 
Hence even when entries have been “‘chopped”’ it is sufficient to study the 
distribution functions associated with 6. 

For purpose of comparison, the known results for ordinary differences 
of orders up to three are summarized below, along with some of the corre- 
sponding distributions for modified divided differences. 

First modified divided difference: 








€1 — €o 
oda Se Sate ie C1 ’ 
(2.28) a a a 
f(wi, t) = (c1/a)? a ip C1 Sts C1 
0 elsewhere. 


If c,; = 1, we have the frequency function for the first ordinary differ- 
ence. 
Second modified divided difference: 


2€0 2e1 2e_; 


We = Mo-1e-1 + Mo,c€0 + Mo,1e1 = — Ge” Gta bd Gita + GS 





Assume Co < ¢; 


(2.29) a [ee (* + ¢; ) fe 





2 2a 
i - a 
4a° » HOS a — ¢1(€o + €3) 
c 
a |- |t\cocx (Co + 1) + aco + 2ac.], 
2a 2a 


,t) = | ae an 
F(won #) . €1(Co + €1) slis Co(Ci + Co) 


(" + 3) (tcoc1)? — 4\t\coc1a + “| 
2 4a’ ; 





, 2a 2a 
if ————~ < [t| < — 
Co(Co + ¢1) — ls Coli 





0 elsewhere. 
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Thus f(we, x) consists of five curves, symmetric about the y-axis, two of 
those curves being linear. If co = c,; = 1, each line shrinks to a point, and 
f (ws, t) consists of three parabolas. 

Third ordinary difference: The following table is taken from Lowan & 
LADERMAN!: 


(iff |e 8a? 

Ald -S 5 +S |. O<l <a 
a. 

> By 4] a |t| < 2a 


(2.30) f(ws,t)=41 ca: af? + ll 
a|s54- 97 +z =|. 2a < |t| < 3a 


1f | , 2a¢ 8a*l 32a ae ye 
a| sat 9 —~9 +97} %#S#S 


L 0 elsewhere. 














There are seven curves in f(ws, t), two of them being straight lines. Let 
us consider the frequency function for the third modified divided difference, 
and write for brevity 


3 
= > Axex; B, = a|A,|. 
k=0 


Making use of (2.27), we can put As; = — (Ao + Ai + Az). If the argu- 
ments x; form an increasing sequence, the sign of A; is independent of the 
magnitude of the intervals (x, — x,_1); hence the sign will be the same as in 
the ordinary third difference, and it is permissible to write 


|As| = |As| — |Ai| + |Adl. 


cos xt[sin (Bot) sin (3B,t) sin (4B2f) sin (4(B. — By + Bo)t) 


G(x, 1) = [BoB:B.(B: — B: + By) 1/16 





By expanding the above into a sum of sines and cosines, it can be verified 
that 


(2.31) G@(x,) =C > {bi (x + D,)* sin [(x + D,)é] 
k=1 


+ by(x — D,)* sin [(x — D,)t]}. 


In (2.31) C is some constant, , =+ 1 if D. + 0, & = 3, if D, = 0, and D, 
assumes the following values: 


(2.32) D; = 0,D2 = Bz — Bi, Ds = Bo, Ds = (Bz — Bi + Bo), 
Ds; = (By = Bo), D, = By, D; = Bo, Ds = Bz 4 Bo. 


Let us consider the special case when 


D,<D2eS Dse DS Ds Del DiS Ds < Bot Bi t+ Bz, 
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and let x be positive. Since 
1 C) 
*f (sin bt/t)dt = Fifb>0, — if b <0, and Oifb =0, 
0 


the terms of (2.31) involving (x + D,)* will contribute, after integration, 
terms that have the same sign for all positive values of x between 0 and 
(Bo + Bi + Bz). On the other hand, the terms involving (x — D,)* will vary 
in sign, depending on whether x is greater than or less than D,. Clearly x can 
fall into any one of the eight regions separated by the inequalities of (2.32), 
and in the most general case when no two D, are equal, the set of terms of 
(2.31) will comprise, after integration, cubic polynomials (or polynomials of 
lower degree) with coefficients that will be different in the several regions. 
Since f(ws, x) is an even function of x, the curve comprising the eight arcs 
will be reflected through the y-axis, with a common central arc. There will 
therefore be 15 arcs to the frequency function for the most general values of 
D,. These shrink to seven arcs in the case of ordinary differences. The labor 
of computing the exact distribution seems to be prohibitive, and alternative 
approximations will be considered. 

It is known that for large enough n, F(w,, t) tends to approach the normal 
distribution with mean 0 and standard deviation o(w,), where 
i 


o(w.) = (= 4s) sia) « o( & 40/12) 


In the above, o(e;) is the standard deviation of e;. It will be instructive to 
examine the probabilities that ws, associated with the third ordinary differ- 
ence, will fall into certain intervals, and to compare them with the prob- 
abilities obtained from the corresponding normal distribution. For ws, 
o(ws) = 1.29099a. The following schedule lists some calculations: 


Theoretical Probability 
Probability Based on Nor- 
Range of ws; Based on (2.30) mal Distribution 
2 o<|ws|< 4a .03687 .04350 
2.40 <|ws|< 4a .00612 .01440 
2.50 <|ws| < 4a .00331 .01236 


It is to be observed that in the third ordinary difference, the normal dis- 
tribution exaggerates the area of the “‘tail’’-end of the distribution. How- 
ever, the discrepancy between the two distributions is no worse than by a 
factor of 2.4 for the first two ranges of the schedule, and agreement is ex- 
pected to be closer in higher differences. 

Assuming that the probability of w; being numerically greater than 2.4¢ 
is approximately the same in higher differences, it might be reasonable to 
tolerate a discrepancy of 2.40. Such a discrepancy is expected to occur once 
in 163 listings of the third difference, according to the exact distribution, and 
if we use the normal distribution as a guide, it may occur once in 69 entries 
in differences of high order. It is costly to examine too many doubtful en- 
tries; often it is much more advantageous to calculate two added decimal 
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places in the entries, so that discrepancies with even higher probability can 
be passed. Much depends on the problem at hand. 

The following schedule lists the value of 2.40 in differences of orders n up 
to 10. Corresponding values of .4>°|A;/a are also included ; for some values of 


‘ . - , P n 
n, .3>.|A;,\a is also tabulated. A, is the binomial coefficient ( | 


n 2.40 AD \A,\a on 2.40 AD \Arla .3>0|Azla 
2 1.697a 1.6a 7 40.588a 51.24 38.4a 
3 3.0984 3.2a 8 78.598a 102.4a 76.84 
4 5.797a 6.4a 9 152.766a 204.84 153.6a 
5 10.9984 12.84 10 297.797a 409.6a 307.2a 


6 21.060a 25.6a 


It should be observed that the simpler function .4>°|A;,|a is close in magni- 
tude to 2.40 if m < 6 and can be used in its place as a basis for a reasonable 
tolerance. In differences of higher order .3>°|A;,|a is close to 2.4¢. 

To what extent can the above schedule be used as a basis for establishing 
a reasonable tolerance for modified divided differences? One way is to com- 
pute a few values of M;,;.4, as a basis for estimating ¢. That may be labori- 
ous. Some qualitative estimates can be obtained from the form of (1.12). 
Let A, denote the coefficient of u, in the ordinary difference of order 2n. 
Then it is clear that 


Mitr Sica oe 
Mo. A; = = i ’ 
P ae (< = Dy il, Cot Cui t ...+ C1 
b sj k 


k # #j 


where s is the greater of k and j and ¢ is the smaller. Thus in a region where 
the values of c;, are consistently greater than unity, the standard deviation is 
lower than that in the ordinary difference, and the tolerance should be more 
stringent than for ordinary differences. On the other hand, if the c; are all 
smaller than unity in a region, then a is larger. Whe-ever the standard devi- 
ation for ordinary differences can be used as a basis, the schedule of 2.4¢ 
offers one means of judging the extent to which an isolated error will be 
masked by smaller errors in neighboring entries. 

Secondary Effects in Round-offs. When ordinary differences are considered, 
all the round-off errors occur in the entires u,, and the process of taking 
differences introduces no other errors. In forming divided differences, how- 
ever, (modified or not) multiplications and divisions are involved, and the 
result is rounded off to a fixed number of decimals. Hence the cumulative 
effect of such roundings must be considered, and we must specify the order 
in which various operations are to be performed. Let g(k,r) denote for 
brevity the rth modified divided difference associated with t,. Then 
g(k, r) is formed as follows: 





g(k,r) = y(k,r)[g(k + 1,7 — 1) — g(k,r — 1)], 
where 
y(k,r) = rw/([tesy — te-y], ifr is even 


y(k,r) = rw/[tesieay) — te-yo—vn], if r is odd. 
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The differences are useful for interpolation or error detection only in the case 
where successive differences fall off in magnitude. Let y(k, 7) be computed to 
the maximum attainable accuracy (depending on the computing machine), 
and then multiplied by g(k + 1,7 — 1) — g(k,r — 1), which in successive 
differences is expected to have fewer significant figures than y(k, r) in the 
useful case. Thus the error due to carrying an inexact y(k, 7) is expected to 
be of a lower order of magnitude than the error in g(k, 7) and we shall 
neglect its consideration. The principal’new error is therefore the rounding 
of g(k, r) to a fixed number of decimals. Thus in each successive difference 
there is introduced a new rounding p;, assumed to satisfy conditions (a) and 
(b). If w, is the error function in the mth difference with exact operations in 
obtaining g(k, r) then the total error due to all types of roundings is 


Pi = Wn -+- Wnr-1 + Wn-2 + see -+ Wi + Pny 


where p; replaces e; in w,_;, 7 > 0. 

It is necessary to have an estimate of V,, as compared with w,. Let us 
consider the case where the standard deviation of w, is used as an estimate, 
and let 

j=€62AP/12: «e = 62. 


Then the standard deviation of V,, is 


1 


ar =. 
w= (Qe + Ea). 


In the special case where the A; are binomial coefficients (i.e., for the 
ordinary difference) it can be verified that 
oy, = Anon 

where d2 = 1.155, ds = ds = 1.183, and for m ranging between 4 and 10, d, 
decreases steadily up to 1.165 for m = 10. Hence if a tolerance of 2.4¢, is 
allowed for ordinary differences, a tolerance of 2.8 ¢, should be reasonable 
for divided differences. If storage space permits, it is of course possible to 
carry two more decimals in the divided differences than in the function, so 
as to lessen the rounding error. If that is done, the added tolerance is not 
necessary. 

The writer acknowledges gratefully the help given by Dr. J. LADERMAN, 
who read the manuscript and offered many valuable suggestions. 

GERTRUDE BLANCH 

NBSINA and 
Consolidated Engineering Corporation 
Pasadena, California 

1 ARNOLD N. Lowan & Jack LADERMAN, “On the distribution of errors in mth tabular 
differences,” Annals of Math. Statistics, v. 10, 1939, p. 360-364. One minor detail of the 


paper is faulty, but the results are correct. In evaluating integrals involving products of 
sines or cosines, it is stated, for example, that 


f di{sin at]/f@"1 = [(— 1)" a/(2n)!] f dt sin at/t. 
_) 0 
This of course is not true, since the left-hand integral does not exist if ” is different from zero. 


However, it turns out that the contribution from the region around the origin in a sum of 
such integrals vanishes, and the result is right. 
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2A. VAN WIJNGAARDEN, Afrondingsfouten MR3, Tevens ZW-(1950)-.001. Math. 
Centrum Rekenfdeling, Amsterdam (in Dutch). See also A. M. Ostrowski. Two Explicit 
Formulae for the Distribution Function of the sums of n Uniformly Distributed Independent 
Variables. Archiv d. Math., v. 3, 1952, p. 3-11. 


RECENT MATHEMATICAL TABLES 


1165[B,F].—H. S. Unter, “On the 16th and 17th Perfect Numbers,” 

Scripta Math., v. 19, 1953, p. 128-131. 

This note contains exact values of 2"~-'(2" — 1) for m = 2203 and 2281, 
numbers of 1327 and 1373 decimal digits. These are the 16th and 17th per- 
fect numbers. Exact values are given also of 2” for n = 560, 2202, 2280 and 
those digits of 2*°* and 2 which are not identical with the corresponding 
digits of the perfect numbers mentioned above. 

The author has informed the reviewer of the fact that the 1023rd digit 
was printed incorrectly: for 32633 read 32638. This substitution occurred 
between page proof and printing and would have gone undetected by any 
author but one having Uhler’s indefatigable perspicacity. 


D. i. o. 


1166[C ].—NBSCL, Tables of 10* (Antilogarithms to the Base 10). NBS 
Applied Math. Series, No. 27, U. S. Gov. Printing Office, Washington, 
1953, viii + 543 p., 19.3 K 26.0 cm. Price $3.50. 

The main table in the work is Table I, a 500 page table of 10* for x = 
0(.00001)1. These 100000 values are given to 10D. The arrangement is in 
four columns of 50 pairs (x, 107) each so that consecutive entries lie one 
under the other making linear interpolation easy. All eleven digits of 10” are 
given in each entry. No differences are given. Linear interpolation gives 9D 
accuracy. The effect of the second difference on the 10th decimal may be 
read from a chart on p. vi. This amounts to at most 7 units in the 10th place. 

Table II is a 15D radix table of 10*. Specifically it gives 10” where 


y=n-10-%, n = 1(1)999, p = 3(3)15. 


From this table 14 figure antilogarithms can be found by multiplying five 
entries together. The table can also be made to serve as a table of common 
logarithms to 15D. 
Table II is similar to that of DEPREz' which gives 13D antilogarithms of 
x = m-10- for 
m = 1(1)999, r = 7(3)13 


in connection with a basic table for 
x = 0(.0001)1. 


Table II will be found very useful in connection with any ordinary radix 
type logarithm table for very precise work. 

Table I is based on Dopson’s? rare table of 1742. The entire Dodson 
table was transferred to punched cards and differenced on a tabulator. 
After correcting errors the table was checked by summing sets of 50 con- 
secutive entries (as a geometric progression). Finally the printed page proof 
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was subjected to additional differencing. No list of errata in Dodson is given. 

Table I is a handy companion to a table of 10 place logarithms since in- 
verse interpolation is avoided in the usual passage from numbers to loga- 
rithms and back to numbers. However since interpolation is only linear (at 
least for 9D work), inverse interpolation presents no greater problem than 
direct interpolation, especially if one has even the smallest type of desk 
calculator. From this there are two alternative inferences: (a) One needs 
no table of antilogarithms; (b) one needs only a table of antilogarithms. 
Perhaps this table will appeal most to those who want to do no interpolation 
whatever. 

This useful volume is the result of work done by the old New York 
Mathematical Tables Project. 

i, T. &. 
1F. DEpPREZ, Tables for Calculating, by Machine, Logarithms to 13 Places of Decimals. 


Berne, 1939. 
2 J. Dopson, Antilogarithmic Canon. London, 1742. 


1167[C ].—NBSCL. Table of Natural Logarithms for Arguments Between 
Zero and Five to Sixteen Decimal Places. NBS Applied Math. Series, No. 
31, U. S. Gov. Printing Office, Washington, 1953, [Reissue of MT 10], 
x + 501 p., 20.1 X 25.9 cm. Price $3.25. 

The original NYMTP Table 10 [v. 3 of the original 4v. ] x = 0(.0001) 5 
is hereby reissued to meet a continued demand. The preface promises also a 
reissue of the fourth volume for x = 5(.0001)10. Although the Introduction 
states that there has been no revision of the tabular content, there has been 
a change in the rule for the indication of the signs of the logarithms. Thus in 
the original edition the logarithm of .0184 is given as 3.995 . . . , the fact 
that this number is actually negative being understood. Now the minus 
sign is printed explicitly. This improvement is carried out with a single ex- 
ception and this occurs at the very first real entry of the table where the 
logarithm of .0001 is given as 9.21034. ... 

One further change may be noted and this refers to the last entry in the 
table. The reader who is familiar with the NYMTP tables will recall that 
arguments are given at the bottom of the page without the corresponding 
functional values. This tantalizing procedure is followed in the present 
volume except at the very end where the editor has relented and has given 


In 5.0000 = 1.6094379124341004. 


Apparently there are no errata known in this monumental table of 1941. 
oo &. 1. 


1168[F ].—E. S. Barnes & H. P. F. SwINNERTON-DyeEr, ‘The inhomogene- 
ous minima of binary quadratic forms,’’ Acta Math., v. 87, 1952, p. 259- 
320. 

Let f(x, y) = ax? + bxy + cy® be an indefinite binary quadratic form 

with real coefficients and discriminant D = b? — 4ac > 0. For any point P: 

(xo, yo), Where xo, yo are real, define M(f;P) to be the lower bound of 


f(x + x0, y+ yo)! 
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taken for all points P and call M(f) the upper bound of M(f;P). Let C be 
the set of points P for which M(f;P) = M(f) and M,2(f) the upper bound of 
M(f;P) over all P not belonging to C. Clearly M2(f) < M(f). If the strict 
inequality holds, M(f) is called an isolated minimum. 

The authors list in their table [p. 315-317] the values of M(f) for forms 
x? — my for all square-free m = 2 or 3 (mod 4), m < 101, except m = 46, 
67, 71, 86, 94 and many corresponding values of M2(f). In many cases com- 
plete sets of incongruent (mod 1) points are given for which M(f;P) = 
M(f) or M(f;P) = M,(f). All minima given in the table are isolated minima. 

A corresponding table is given for forms f = x* + xy — }(m — 1)y 
where m = 1 (mod 4), m is square-free and not greater than 101, except 
m = 57 and 73. 

Sources for the results listed, many in the accompanying paper, are given. 


B. W. JONEs 


Univ. of Colorado 
Boulder, Colo. 


1169[ F ].—A. GLopEN, Table des Solutions de la congruence x‘+1=0 
(mod p) pour 600000 < p < 800000. Luxembourg, 1952, published by 
the author, rue Jean Jaurés, 11, Luxembourg, 22 p., 29.8 K 21.0 cm., 
mimeographed. Price 120 francs belges. 

This table is identical with that described as UMT 158, MTAC v. 7, p. 

108 except that the appended factorizations there referred to are wanting. 


1170[F,L ].—A. vaAN WIJNGAARDEN. “On the coefficients of the modular 
invariant J(r),’’ K. Ned. Akad. v. Wetensch. Proc., s.A, v. 56, 1953, p. 
389-400. 

KLEIN’s fundamental modular invariant J(7r) has an expansion in the 
form 


12°J(r) = x"! + 744 + 196884x + 21493760x? +--- 


x 


> c(n)x" 


n=—1 


in which x = qg = e’*'". These coefficients are positive integers whose prop- 
erties have only recently been investigated. The present paper gives a two 
page table of c(m) for m = — 1(1)100. These values are quite large, c(100) 
having 53 digits, and the table represents a relatively large amount of com- 
puting. Previous tables are those of Berwick! for n < 7 and ZUCKERMAN? 
for n < 24. 

The intimate connections between J(7) and other elliptic functions pro- 
vide a variety of methods for the computation of c(m). That none of these 
methods can be too easy is pointed out by the author who remarks that “the 
coefficients grow very rapidly with m and the digits have to come from some- 
where.”’ 

Zuckerman exploited a connection between c(m) and the number of parti- 
tions of 25”. This becomes ineffective for m = 25 since Gupta’s* tables of the 
partition function extend only to n = 600. 
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LEHMER*‘ had proposed the formula 


n 


Dd c(k)r(n — k) = 720{91011(m) + 600r(m)} /691 


k=-1 


and VAN DER Po -' the more elegant formula 
(1) > ke(k)r(n — k) = 24043(n). 
k=—1 


These formulas are based upon the connection between J(r) and the 
Weierstrassian discriminant 


i) 


Dd r(n)x" = x II (1 — xn)™ 


n=1 n=1 


whose coefficients, known as Ramanujan’s function, are now tabulated® to 
n = 2500. Formula (1) was used by the author up to m = 50. At this point 
another more elaborate method based on the relation 


27] (7) = 2 (08 + 638 + 648) (65-8 + @3-§ + 04-8) 


between J(7r) and Jacobi’s theta functions, was used to recompute and ex- 
tend the table to m > 100. Various congruence properties of c(m), including 
an interesting new one modulo 71, were used to check the calculations. 

The table should be a valuable tool for further research on Klein’s in- 
variant. 

2. . §. 

1W. E. H. Berwick, “An invariant modular equation of the fifth order,”’ Quart. Jn. of 
Math., v. 47, 1916, p. 94-103. 

2H. S. ZuCKERMAN, “The computation of the smaller coefficients of J(r),’’ Amer. 
Math. Soc., Bull., v. 45, 1939, p. 917-919. 

3H. Gupta, ‘A table of partitions, I, II,’ London Math. Soc., Proc., v. 39, 1935, p. 
142-149, v. 42, 1937, p. 546-549. 

4D. H. LeuMer, “Properties of the coefficients of the modular invariant J(r),’’ Amer. 
Jn. of Math., v. 64, 1942, p. 488-502. 

5B. VAN DER Pot, “On a non-linear partial differential equation satisfied by the loga- 
rithm of the Jacobian theta-functions, with arithmetical applications, I, II,’”’ K. Ned. 
Akad. v. Wetensch., Proc., s.A, v. 54, 1951, p. 261-284. 

6 See MTAC, v. 4, 1950, p. 162, UMT 101. 


1171[1].—H. E. Sauzer, Tables of Coefficients for the Numerical Calculation 
of Laplace Transforms. NBS Applied Math. Series, No. 30, U. S. Gov. 
Printing Office, Washington, 1953, ii + 36 p., 20.1 & 25.9 cm. Price 25 
cents. 

These tables are intended to be used in the approximation of the trans- 
form 


F(p) - f e-P'f (t)dt 


by means of the sum 


n—1 


F(p) = X A;f(i) 
i=1 
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where 
A; = A;™ (p) 


depend on both # and , the number of ordinates. The main tables are ar- 
ranged by m which takes the values 2(1)11. The values of p are .1(.1)” — 1 
up through » = 7. For m = 8 and 9 the interval is .2 and for m = 10 and 11 
it is 1. Values of the A’s are given to 9S. There are auxiliary tables as follows. 

Table II (p. 27—36) gives values of n!/p"*! for n = 0(1)10, p = .1(.1)10 
to 8s. This function is the Laplace transform of ¢” and the table is intended 
for use in transforming polynomials up through those of degree 10. 

The function (m — 1)!p"A;™(p) is a polynomial in p with integer co- 
efficients. These polynomials are listed in the introduction (p. 7-8) together 
with the corresponding Lagrange interpolation coefficients (polynomials in f) 
for m = 2(1)11. 

Four explicit examples are worked out including one in which the error 
is estimated. 

These tables should be quite useful in dealing with functions whose 
Laplace transform is unfamiliar or intractable. 


D. H. L. 


1172[K ].—Hrirojiro Aoyama, ‘On a test in paired comparisons,” Tokyo 
Inst. Math., Annals, v. 4, 1953, p. 83-87. 
Let each of ” persons, no two of which have the same occupation, rate his 
own occupation in comparison with each of the others. This gives rise tok = 


2 

those cases in which each subject in a pair rates his own occupation higher 
than the other, —1 in those cases in which each subject rates his own occupa- 
tion lower than the other, 0 in all other cases. Let S be the total of the & 
scores. This is proposed as a test criterion for the null hypothesis that each 
subject is as likely as not to rate his own occupation higher than any other 
of the ” occupations. Under this model for k throws of a pair of coins, S is 
the excess of the number of cases of both heads over the number of cases of 
both tails. A table_is given for Pr(|.S| 2 So) for m = 3(1)9 to at least 4D 
always giving at least 2S. 


n ‘ . . ‘ , 6 
( paired comparisons, to which we assign the following scores: +1 in 


c eS 


1173_K ].—joseru Berkson, “‘A statistically precise and relatively simple 
method of estimating the bio-assay with quantal response, based on the 
logistic function,’’ Amer. Stat. Assn., Ju., v. 48, 1953, p. 565-599. 

This paper sets forth very clearly and succintly the author’s ‘“‘minimum 
logit chi-square’’ method in bio-assay. To facilitate the necessary computa- 
tions three tables are provided. Table 1 gives / = In [p/(1 — p)] to 5D for 
p = .001 (.001) .999. Table 2 inverts Table 1, giving p to 5D for / = 0(.01)- 
4.99. Table 3 gives the weights needed in the logistic calculation, namely 
p(1 — p) and p(1 — p)l, to 4D, for p = 0(.001)1. 

J. L. Hopces, Jr. 

University of California 

Berkeley, California 
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1174(K ].—K. N. CHANDLER, ‘‘The distribution and frequency of record 
values,’”’ Roy. Stat. Soc., Jn., s.B, v. 14, 1952, p. 220-228. 


The frequency and interval of occurrence of record values is of con- 
siderable interest whether we are dealing with weather data or sampling 
other types of populations. The lowest record value at any given time is de- 
fined as that member of the sample which is less than or equal to all previous 
members and the greatest record value similarly is defined as that value 
which is greater than or equal to all previous values. In this paper Chandler 
considers the random series x,, u = i, 2, 3, etc. with the first record value X, 
equal to x;. Let X; be the first occurring value of x which is less than X;_,, 
4 = 2,3, etc. Then X, is called the rth lower record value (similar considera- 
tions apply to the higher record values). Assuming that the distribution 
function of x is either normal or rectangular, Chandler derives the probabil- 
ity distribution of the rth lower record value X, (or rth higher record value), 
the distribution of the serial number, u,, of the lower record value and also 
the probability distribution for the interval of occurrence (number of ob- 
servations) between the 7th lower record value and the (r — 1)st lower 
record value. 

Table 1 and Table 2 of the paper give the .005, .01, .1, .5, .9, .99 and .995 
probability points for the distribution of X, for the normal distribution in 
standard units to 3D and for the rectangular distribution to 4S for r = 
2(1)9 in both cases. Table 3 gives the probability that the serial number, 
u,, of the lower record value, will be < to 6D for r = 3(1)9 and m = 3(1) 
30(5) 60(10) 100, 200, 500, 1000, 2000, 5000, 10000, 20000, 50000, 100000, 
200000, 500000, 1000000. Table IV gives the probability that the interval 
of occurrence, u, — u,—;, the number of observations between the 7th lower 
record value and the (r — 1)st lower record value 2 » to 6D for r = 2(1)9 
and the same set of values of m as in Table 3. The serial numbers u,, and their 
differences, u, — u,_;, do not have finite means and are independent of the 
parent population provided it has finite or zero probability density at all 
points! 

F. E. GRUBBS 
Ballistic Research Laboratories 
Aberdeen Proving Ground, Md. 


1175[K ].—W. J. Dixon, ‘Power functions of the sign test and power 
efficiency for normal alternatives,’ Annals Math. Stat., v. 24, 1953, p. 
467-473. 

The two sided sign test has long been used as a non-parametric test of 
significance. In using such a test one generally compares the number of 
changes in sign that occur in his observations with the expected number 
under a given hypothesis. In using such a non-parametric technique, one is 
interested in its power as well as how this power compares with the power of 
corresponding parametric tests. This paper contains tables which help the 
research worker answer such questions. 

Tables I and II give respectively the power of the sign test to 5D for the 
5% and 1% level of significance for sample sizes N = 5(1)20(5)50(10)100 
in testing the null hypothesis that the proportion of objects in the population 
is .5 against varying proportions, p = .05(.05).95 in the alternative popula- 
tion. 
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To compare the power of the sign test with normal alternatives the 
author introduces a power efficiency function which gives the power effi- 
ciency of the sign test as compared with the normally based test for each 
alternative. Tables for selected values of the parameters and representative 
curves of this power efficiency function are given. The results exhibited by 
this function are compared with various asymptotic and approximate 
efficiency estimates that have been obtained by various authors in the past. 


C. F. Kossack 
Purdue University 
Lafayette, Indiana 


1176[K ].— Benjamin Epstein & Mitton Soset, “Life testing,” Amer. 
Stat. Assn., Jn., v. 48, 1953, p. 486-502. 


For a characteristic X with density function 


f (x30) = (exp (— x/6))/0 


the maximum likelihood estimate of @ based on the first r order statistics of a 
sample of 1 is 


- = [x10 +: ’ + Xrn + (n nay 1)X-r,n]/P. (r 4 n) 


Furthermore, 274,,,,/@ is distributed as x? with 2r degrees of freedom. Since 
this result is independent of nm, an appreciable saving of time can be made in 
situations, such as life testing, where the observations become available in 
order by basing the estimate of 6 on the first r available results from a larger 
sample. Table I tabulates to 2D the ratio E(X,,.)/E(X,,,) of the expected 
time to obtain the first r results from a sample of ” to the expected time to 
obtain all r results from a sample of r for r = 1(1)5, 10 and m = 1(1)5(5)20. 
A simple derivation of previously known forms for E(X,,,) and Var(X,,.) is 
given. 

It is also shown that the region of rejection for the best test of H,(@) = 6; 
against the alternative H2(@) = 02 < @, is given by 6, < C. Table II 
tabulates for 6;/@2 = 1.5(.5)3(1)5(5)10; a, 8 = .01 and .05 the minimum 
value of r, and the corresponding upper and lower limits for C/@, to 4D, such 
that the errors of Type I and Type II will be less than a and 8, respectively. 
This table represents a rearrangement and, in some respects, an extension of 
tables published by E1rsENHART! which (in the notation of the present paper) 
gives for fixed 7, a, and 8 the maximum value of 6;/@2 such that the test con- 
ditions are satisfied. It is also directly related to other tables dealing with 
the operating characteristics of the one-sided tests of a hypothetical variance 
oo” against the alternatives o? < a9. 

C. A. BENNETT 
General Electric Company 
Hanford Atomic Products Operation 
Richland, Washington 


1C. E1sennart, M. W. Hastay, W. A. WALLIs, editors, Selected Techniques of Statistical 
Analysis. New York and London, 1947, ch. 8, p. 267-318. 


1177[K ].—H. L. Jones, “Approximating the mode from weighted sample 
values,”” Amer. Stat. Assn., Jn., v. 48, 1953, p. 113-127. 
This paper presents a method of estimating the population mode by a 
weighted sum of the order statistics of a sample. The derivation of the values 
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of the weights is based on an approximation to the maximum likelihood solu- 
tion. The parametric form of the population is assumed to be known. The 
method is of a general nature but does not necessarily yield a reasonable ap- 
proximation to the mode unless certain favorable conditions are satisfied. 
The case of a #-distribution with varying degrees of kurtosis is analyzed to 
illustrate application of the method. Table 1 contains weights to 2D to be 
used in approximating the mode of a ¢-distribution for m (sample size) = 
3(1)10 and ag = ps/p2? = 3(.5)5, 6,9, @. 

J. E. WALsH 
U. S. Naval Ordnance Test Station, Inyokern 
China Lake, California 


1178[K].—R. Kamar, “On the mean successive difference and its ratio to 
the root mean square,”’ Biometrika, v. 40, 1953, p. 116-127. 
Given a sequence of m normal variates {x;} with common mean and 
variance. The mean square successive difference is 


n—) 
d = (n —1)7*D [Xi — Xigil. 
i=1 
Table 1 presents the standard deviation, 8; and 2 of d/c, for m = 3(1) 
10(5) 30, 40, 50 to 4D. Table 2 presents approximate upper and lower .5, 1, 
2.5 and 5% percentage points of d/o, to 2D, using a Pearson Type I curve; 
exact results are given for » = 3. 

The ratio W = d/s, where s is the sample standard deviation, is also 
considered. Table 4 presents the mean to 4D, the standard deviation to 4D, 
8; to 3D and 6: to 2D of W for nm = 5(5) 30, 40, 50. Upper and lower 0.5, 1, 
2.5 and 5% points for W are given in Table 5 for m = 10(5)30(10)50 to 2D. 

R. L. ANDERSON 
North Carolina State College 
Raleigh, North Carolina 


1179[K].—Tosio KitaGAwA, TEISUKE KITAHARA, YUKIO Nomacur & 

Nosuo WATANABE, ‘‘On the determination of sample size from the two 

sample theoretical formulation,’’ Bulletin Math. Stat., v. 5, 1953, p. 

35-46. 

The authors consider a two-sample procedure alternative to that of 
Stein!” for determining a confidence interval of fixed length 2d for the mean 
of a normal population with unknown variance. They give a table for the 
following function connected with the probability distribution of the second 
sample size: 

b(n) 


I(n2;m|d?o~,a,8) = f 


b(ne—1 


e(ng) 
dnlsizoddss f $n, (S2;o)ds2, 
where 
b(m) = dn*{tp1(a)}-*{ Fai (8)}-, 
c(nm) = dn*{t,1(a)}—, 


tn-1(a) is the 100a-percentage point of the é-distribution with n — 1 degrees 
of freedom, Fri (8) is the 1006-percentage point of the F-distribution with 
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n — 1 and m, —1 degrees of freedom, 


ni ; 
= (n; —1)7* D (wy — %)??, * = 1,2, 
j=1 
is the square root of the unbiased estimate of the population variance given 
by the first, respectively, second sample, of sizes m; and m2, and ¢n;,(s;30) 
denotes the density function of s;. The tables are for only the selected com- 
binations of the arguments m,; = 10, 15, 21, 25, 31; mz = 3(1)60; a = .01, 
.05; 8 = .01, .05; da? = .75, .5, .25. 
Jutius LIEBLEIN 

NBSSEL 


1W. G. Cocuran, Sampling Techniques. New York, 1953, p. 59-60. 
2B. M. SEELBINDER, “On Stein’s two-stage sampling scheme,” Ann. Math. Stat., v. 24, 
1953, p. 640-649. 


1180[K ].—J. Lerévre, ‘Application de la théorie collective du risque a la 
réassurance ‘Excess-Loss,’”’ Skandinavisk Aktuarietidskrift, v. 35, 1952, 
p. 161-187. 


This paper contains a table of values of 


B,(u) = [~ tedgn(e 


in which ¢(¢) is the cumulative normal frequency function in standard units. 
Values are given for n = 0, 3, 4, 6 for u = 0(.1)3 to 5D. 
a a See 


1181[K].—Jack Mosuman, “Critical values of the log-normal distribu- 
tion,”” Amer. Stat. Assn., Jn., v. 48, 1953, p. 600-609. 


The three-parameter log-normal distribution function may be written as 


(2x)-* 1 \ x— *) 

c(x — a)? — 35(tos b } 

In terms of parameters a, b, and c, the mean, variance and skewness (third 
standard moment) of this distribution may be expressed as un = bw! + a, 
o? = Bw(w — 1), and as = + (w — 1)*(w + 2) where w = exp ¢ and az 
takes the same sign as 6. The author tabulates selected critical values rg of 
the standardized log-normal variate such that P(r = rs) = 8, where r = 
(x — u)/o. Tabulations are to 3D for 6 = .005, .01, .025, .05, .10, .90, .95, 
.975, .99, and .995 with as = 0(.05)3. Accuracy to within one or two digits 
in the last decimal is claimed for all table entries and the author indicates 
that three point Lagrangian interpolation will give similar accuracy for 
intermediate values of as. 





f(x) = 


A. C. CoHEN, JR. 
University of Georgia . 
Athens, Georgia 
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1182[K ].—K. R. Narr, “Tables of percentage points of the ‘Studentized’ 
extreme deviate from the sample mean,” Biometrika, v. 39, 1952, p. 
189-191. 


Let x,(v = 1(1)”) be the rth ordered variate in a sample of size m taken 
from a normal population with unknown standard deviation o. If an estimate 
s, of o is available with r degrees of freedom, independent of the sample, the 
author suggested the use of the Studentized extreme deviation (x, — )/s, 
or ( — x,)/s, as a test criterion for a single outlier. In a previous paper! he 
gave the lower and upper 5% and 1% points of this deviate. These tables 
have now been extended to cover four more percent points, namely 10, 2.5, 
.5 and .1%. Table 1A gives the 6 lower percentage points mentioned to 2D 
for m = 3(1)9 andr = 10, 15, 30, ©. Table 1B gives the same six upper per- 
centage points to 2D for m = 3(1)9 andr = 10(1)20, 24, 30, 40, 60, 120, ~. 

E. J. GUMBEL 
Columbia University 
New York, New York 


1K. R. Narr, ‘The distribution of the extreme deviate from the sample mean and its 
Studentized form,” Biometrika, v. 35, 1948, p. 118-144. 


1183[K ].—E. S. Pearson & H. O. HARTLEY, “Charts of the power function 
for analysis of variance tests, derived from the non-central F-distribu- 
tion,’’ Biometrika, v. 38, 1951, p. 112-130. 


Let u;(¢ = 1, 2,---, v) be » normally distributed independent variables 
with unit variance and zero mean and let a;(i = 1, 2,---, v) be » fixed con- 
stants, then the distribution of 


x” _ > (u; + a;)’, 


i=1 
which is a Bessel function, is called the non-central chi-square distribution 
with v degrees of freedom and \ = 5 a? is called the non-centrality param- 


1 
eter. If x1” is such a value with »; degrees of freedom and x;* another inde- 
pendent central chi-square with v2, then F’ = (x1v2)/(x2'v1), called the non- 
central variance ratio, has a known distribution. Its numerical values are 
obtained with the help of the incomplete 8 functions. The probability 6(A|a, 
v1, ¥2) = Pr (F’ > F,) regarded as a function of ) is the power function of 
the analysis of variance test with significance level a. On the basis, mainly, 
of TaNG’s table! eight charts are given corresponding respectively to v1 = 
1(1)8 for B at the levels a = .05 and a = .01. Here the non-centrality param- 
eter ¢ = (A/(v; + 1))? is used as the abscissa instead of \ on a linear scale. 
Each chart gives two families of eleven power curves corresponding to vz = 
6(1)10, 12, 15, 20, 30, 60, © for the two values of a. The use of a logarithmic 
scale for the ordinate 8 straightens the curves and expands them in the re- 
gion of high power .80 S 8 < .99. The B grid, .1(.1).5(.05).7(.02).90(.01).99, 
and the ¢ grid of .2 for a = .01 and of .1 for a = .05 are sufficiently fine to 
allow interpolation by sight. The calculation of g is shown for the one way 
classification into k groups with m observations in each, for double classifica- 
tion with one observation in each cell and for the latin square arrangement. 
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The use of the charts is explained for the analysis of the effect of machine 
variations on the standard deviation of manufactured bulk product and of 
the effect of personal factors introduced in routine tests. 

E. J. GuMBEL 


Columbia University 
New York, New York 


1P. C. TanG, “The power function of the analysis of variance test with tables and 
illustrations of their use,” Stat. Res. Memoirs, v. 2, 1938, p. 126-157. 


1184[K ].—Frank Proscuan, “Confidence and tolerance intervals for the 
normal distribution,” Amer. Stat. Assn., Jn., v. 48, 1953, p. 550-564. 


The author presents an excellent summary of confidence and tolerance 
intervals for the normal distribution for the various combinations of known 
and unknown mean and standard deviation. Let x be normally distributed 
with mean yz and standard deviation oc. Define 


LQaSX/e, 2 eT (ke: - BY/e- pd. 


i=1 t=1 


Let m represent either u or X; let s.d. represent either o or s. Then either 
confidence interval statements or tolerance interval statements may be 
made about m + k s.d., where the value of k depends on the particular type 
of interval and whether or not » and o are known or unknown. All tables of 
k;, i = 1(1)9 are given for m = 2 (1) 30, 40, 60, 120, ~, to 3D. 

Given o known, the 50% confidence interval for pu is given by X + hk; a, 
k, = .6745//n; if both uw and o¢ are unknown, the 50% confidence interval is 
given by X + ks; if both w and o are unknown, X; + kss provide 50% con- 
fidence interval for the second sample mean X2, when 1 = m2. 

The next k;,z = 4(1)7, refer to tolerance limits. If both » and o are known 
u + kyo, kg = .6745 provide tolerance (probability) limits such that the 
proportion p of the population included by the interval is .50. In case » and 
o are unknown the factors ks,, in the tolerance limits X + ks,.s are given for 
a = .50, .75, .95, .99, .999, i.e. the average p contained in X + ks,.5 will be a. 
If u is unknown and o known, tolerance limits are found by use of X + keo, 
and k, is tabulated for a = .50; for » known, o unknown, tolerance limits 
are given by uw + ks and k; is tabulated for a = .50. 

The final cases kg and ky refer to confidence statements about tolerance 
limits. BowKER! has tabulated the values of k such that the probability is 7 
that X + ks will include p or more of the population. In Bowker’s tables p 
and y are listed for all combinations of .75, .90, .95, .99 and .999, » and o 
unknown. The author assumes » known and o unknown and gives values of 
ks for p = y = .50; and uw unknown, o known, ky for p = y = .50. 

L. A. AROIAN 
Hughes Research and Development Laboratories 
Culver City, California 


1C. E1rsENHART, M. W. Hastay, W. A. WALLIS, editors, Selected Techniques of Statistical 
Analysis. New York and London, 1947, ch. 2 by A. H. Bowker, p. 102-107. 
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1185[K ].—S. Rusuton, “‘On sequential tests of the equality of variances of 

two normal populations with known means,” Sankhyd, v. 12, 1952, 

p. 63-78. 

Tables are given for sequential tests of the equality of variances in two 
samples from normal populations using sums of squares (a) about population 
means, (b) about sample means and using ranges of four and eight observa- 
tions. U is the ratio of sums of squares from the two samples. Sampling is 
continued unless U < U’ (accept o1 = o2) or U > U” (accept o1 = 602). 
Tables list U’ and U”’ for degrees of freedom = 1(1)10(2)20,25,30 and 6 = 
1.5, 2.0, 3.0 for all combinations of a and 6 = .01, .05, .10. The test using 
ranges is based on R,, the ratio of sums of ranges from the two samples. The 
tables list R,’ and R,”” for the test as above and for the same range of param- 
eters except that a and @ are restricted to .01 and .05 and m refers to the 
number of ranges of 4 or 8 observations so that the sample sizes are 4n and 
8n. 

W. J. Dixon 
University of Oregon 
Eugene, Oregon 


1186[K].—E. S. Smitru, Binomial, Normal and Poisson Probabilities. 
Published by the author, Box 224C, RD2, Bel Air, Md., 1953. 71 p., 
21.6 X 27.9 cm. $2.50. 

This is a small set of tables and charts, centering on the cumulative 
binomial distribution, with extended discussion of actual (some novel) as 
well as alternative methods of computation. Many of the tables shown are 
obtainable from larger tables constructed by others (referred to in the text). 
Tables are: 


(1) Maximum Poisson probabilities »;(x,a) for various a (x, of course, 
= 0, a, a — 1), a monotone function falling from .999001 (at x = 0) for 
a = .001 to .039861 for a = 100 to 6D 


(2) n! and log n!, n = 0(1)200 10S and 10D 
respectively 

(3) log n, = 1(.01)10 10D 

(4) C,” and log C,", n = 1(1)50 5S and 5D 
respectively 

(5) e-, = 0(.001)1(1)100 10D 

(6) Bi(c,n,p), nm = 1(1)20; p = .01(.01)5 5D 


where B(c,n,p) bad > C."p*(1 — p)"* 


(7) Normal (Gaussian) integrals [‘sora, density functions ¢(¢), second 
derivatives ¢® (¢), ¢ = 0(.01)4, to 5D. 

(8) p(c,a), c=1(1)22; a=.001(.001).01(.01).1(.1)2(1)10 and 
a = 10(10)100 to 5D, 


c 
= = .1(.1)2.2; 


where p(c,a) = = e~*a*/x! 





1] 
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The cumulative binomial probabilities are obtained (a) directly, (b) from the 
Poisson cumulatives (singly or the two-term Gram-Charlier Type B), (c) 
from the normal (singly or the two-term Gram-Charlier Type A, with or 
without remainder modifications) depending on the ranges of the binomial 
parameters. There are a number of useful charts, not, to my knowledge, to be 
found elsewhere. These include (1) the probability of 2 successes in 7 trials, 
with constant probability of success in single trial, (2) the expected number 
(np) of binomial successes, (3) (most useful of all) a set of charts, some 
expressed in terms of correction factors, comparing the normal, binomial, and 
Poisson probabilities for various ranges of the parameters. 

H. A. FREEMAN 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 


1187(K ].—H. Tuer, “On the time shape of economic microvariables and 
the Munich business test,’’ Inst. International Stat., Revue, v. 20, 1952, 
p. 105-120. 


This paper contains a table giving the frequency distributions of the 
difference between two independent random variables each obeying Poisson 
distributions with means m and u respectively. Values are given to 3D for 

1 
mw = 2 R, 2, 4. 


C=. 


1188(K ].—H. Uranist, “The distribution of statistics drawn from the Gram- 
Charlier Type A population,” Bull. Math. Stat., v. 4, 1950, p. 1-14. 


From the expansion of a frequency function in Gram-Charlier Type A 
series, the author obtains early terms of a similar expansion of the distribu- 
tion of the é-statistic for a sample of m. For both the one-tailed and two- 
tailed cases, four coefficients of this series are tabled to 6D for = 5, 10, 15, 
21 for values of the argument u = (1 + #/[m — 1])-! = .05, .1(.1).9, .95, 1. 
With the aid of these tables it is possible to estimate the tail probabilities 
for given values of 81, B2, 83, and 84, and this is done for several examples. 
The results are not entirely comparable with those of GAYEN,' in that the 
latter employed the Edgeworth series. 

J. L. Hopces, Jr. 
University of California 
Berkeley, California 


1 A. K. GayeEn, “The distribution of ‘Student’s’ ¢ in random samples of any size drawn 
from non-normal universes,” Biometrika, v. 36, 1949, p. 353-369. 


1189[K ].—J. WeESTENBERG, ‘“‘A tabulation of the median test with com- 
ments and corrections to previous papers,’’ K. Ned. Akad. v. Wetensch., 
Proc., v. 55, s.A, 1952, p. 10-15. 


Let Xi,---, Xw, and Y;,---, Yy, be samples of Ni, and Nz observations 
drawn at random from populations having continuous distribution functions 
F(x) and F.(y) respectively. Let 5 be the number of observations belonging 
to one of the samples that lies between the median of that sample and the 
median of the combined sample. The hypothesis Hp that F; = F, is rejected 
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when 6 exceeds a critical value 59. Tables of 59 to 1D are given for significance 
levels .001, .005, .01, (.01), .05 when considering two-sided alternatives to 
Hy and .0005, .0025, .005, .01, .015, .02, .025 when considering one-sided 
alternatives to Ho for Ni, N2 = 6, 10, 20, 50, 100, 200, 500, 1000, 2000. 


Cyrus DERMAN 
Columbia University 


New York, New York 


1190[K ].—J. W. WuitFIELD, “The distribution of total rank value for one 
particular object in m rankings of m objects,’ British Jn. of Stat. Psy- 
chology, v. 6, 1953, p. 35-40. 

The problem considered is, essentially, that of the distribution of the 
sum of m independent random variables, each with uniform distribution on 
the integers 1, 2,---,. The distributions, of course, are symmetric about 
3m(n + 1). Three pages of tables give the lower halves of the cumulative 
distribution functions to 5D for m = 2 (1) 8 and m = 3 (1) 8. 


LEo Katz 
Michigan State College 
East Lansing, Michigan 


1191[L].—C. Doms, “Tables of functions occurring in the diffraction of 
electromagnetic waves by the earth,”’ Advances in Physics, v. 2, 1953, 
p. 96-106. 


The tables given in the paper are related to the Airy integral 


Ai(z) = x f cos (4/ + szt)dt 
0 


in the complex plane. 
The numbers — a, being the zeros of Ai(z), the author puts 


i[— a, + y exp (xi/3)] 


‘ A 
fa(y) = exp (An + ttn) = Ai’ (— ay) exp (xi/3) “ 





. the numbers 5, being roots of the equation 
Ai(z) = tAi’(z) exp (— 52i/12), 
in which ¢ is a parameter, he also puts 


£, + inn = 6, exp (xi /6) 
and 


exp (Yn + 5,) = 347! exp (— Sai/12)[Ai’ (b,) 
xX [1 — #b, exp (— 5xi/12)T-. 


Tables of the Airy integral for real argument have been reviewed in 
MTAC, v. 2, p. 302-305, RMT 413, and for complex argument, in v. 2, p. 
309, RMT 420. 

The author tabulated the functions f,(y) for » = 1(1)5 in 1942, under 
the guidance of J. C. P. MiLter. These tables were subsequently checked, 
corrected, and sub-tabulated by the Mathematics Division of the National 
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Physical Laboratory of Great Britain. Photostatic copies are available from 
H. M. Nautical Almanac Office, Great Britain. The computations of 
Admiralty Computing Service are described in MTAC, v. 2, p. 35, RMT 260. 

Table 1 of the present paper gives 3D values of A, and yu, for m = 1(1)5, 
y = 0(.2)3(1)10 and 3D values of A, and u, + Fy! for m = 1(1)5 and y = 
10(10)100. 

Table 2 gives 3D values of & and », for m = 1(1)5, ¢ = 0(.1)1, f° 
= 1(—.1)0. 

Table 3 gives 3D values of 7, and 6, for m = 1(1)5 and for ¢ ranging 
from 0 to  ; the number of selected values of ¢ varies with n. 

There is a brief description of the computations, and an indication of the 
application of the functions tabulated here. 

A. E. 


1192[L].—M. Masuixo, Tables of Generalized Exponential-, Sine-, and 
Cosine Integrals Ei(x + iy), Si(x + iy), Ci(x + ty). Numerical Compu- 
tation Bureau, Tokyo, Japan, Report No. 7, March 1953, 43 p. 

Let 


I(z) = f te—‘dt. 


For z = te* (O< ~< 5) put 
I(ge*) = C.(€) — +S.(€). 
For z = ~ee(n 0.2) put 


I (: e) = — a(n) exp (7.(n)). 


The report contains values of 


(a) C.(€) + log — and S.(£) to six decimal places with second differ- 
ences for — = 0(.05)5.00, a = 0°(2°)60°(1°)90°. 

(b) Aa(m) to six decimal places, ,(y) to five decimal places, each with 
second differences, for the above range of a and 9 = 0(.01).20. In defining 
the ranges, the reversals of the inequality signs are presumably misprints. 


The present table is therefore concerned with the exponential integral 
for complex argument, not the generalized exponential-integral which has 
been tabulated by the Harvard University Computation Laboratory 
[LMTAC, v. 4, p. 92-93]. It breaks new ground in covering, for complex 
arguments in polar form, the whole first quadrant. 

No statement is made as to the accuracy of the table. Spotchecking a 
few values against the forthcoming NBS Tables of Exponential Integrals for 
Complex Arguments (cartesian form) did not reveal any discrepancy. 

Everett’s interpolation formula involving second differences gives the 
maximum attainable accuracy in both directions. In order to facilitate inter- 
polation, a four decimal place table of Everett second-order interpolation 
coefficients for arguments at intervals of .01 is also given. 

I. A. STEGUN 
NBSCL 
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1193[L].—NBSCL, ‘“‘Struve function of order three-halves,’’ NBS, Jn. cos 
Research, v. 50, 1953, p. 21-29. ing 
Tables of 
Qn \t Mic 
hy(x) = (2 Hi, (x) Ka 
2 2 cos x 
di & iit ° 11¢ 
etre 2 ( sins + = ). 


to 10D, with second central differences, modified when necessary (modifica- 

tion being indicated by a letter C placed after the entry); x = 0(.02)15. 

“The values are expected to be correct to within one unit of the last place.” 
The entire computation of this table was done on the SEAC, under the 

supervision of ETHEL MARDEN, by KATHRYN CHRISTOPH, ANNE FUTTER- 

MAN, RENEE JASPER, SALLY TSINGOU, and BERNARD URBAN. for 
The table is also available on IBM cards. 


A. E. I, 


1194[L ].—Hersert E. Sauzer, Ruta Zucker, & RutH CAPuANo, “Table 
of the zeros and weight factors of the first twenty Hermite polynomials,” | 11 
Jn. Research, NBS, v. 48, 1952, p. 111-116. 
The definition of Hermite polynomials used in this paper is 


see 


‘ 





oe 
= o>. 


H,(x) = (— 1)re* 3 


x; is the i-th positive zero of H,(x), 

a = h2mn CH! (xs) fo 

is the corresponding Christoffel number, and . 
BO = a exp [(xi)?]. 


The present paper gives 15D values of x; and 13S values of a; and 
Bs for nm = 1(1)20, ¢ = 1(1)m. A list of references (28 items) is appended. 
Other tables of zeros of Hermite polynomials are referred to in MTAC, 

v. 1, p. 152, RMT 131; v. 3, p. 26, RMT 466; v. 3, p. 416, RMT 619; v. 3, 
p. 473, RMT 641; v. 6, p. 232, RMT 1034. The present tables were compared 
with those of Re1z [RMT 466], GREENwoop & MILLER [RMT 619] and 
Korat [RMT 641], and with the Harvarp tables [RMT 1034]. ‘ 
A. E. 


1195[L ].—K. M. SrecELt, J. W. Crispin, R. E. KLEmnMAN & H. E. HUNTER, 
“Note on the zeros of (dP,,, (x)/dx)|z_2," Jn. Math. Phys., v. 32, 1953, 
p. 193-196. 
The authors apply the identical technique of a previous paper to obtain 
approximate values m; such that dP™,,;(xo)/dx = 0 and values of 





[=P as) de . 


[see MTAC, v. 7, p. 183]. Again, the theory is demonstrated for xo= 
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cos 165°. For this case, the first 15 approximate values of m; and correspond- 
ing values of the integral are tabulated. 


YuDELL L. LuKE 
Midwest Research Institute 
Kansas City, Missouri 


1196[L].—N. B. SLater, “Gaseous unimolecular reactions: theory of the 
effects of pressure and of vibrational degeneracy,” Roy. Soc. London, 
Philos. Trans., s.A, v. 246, 1953, p. 57-80. 


Table 3 (p. 71) gives values to 2 or 3S of 
I,,(0) = Crim +i f x™e-*(1 + O-1n™)—“Idx m = (n + 1)/2 
0 


for n = 3(2)13 and some of the values logio? = — 2(1)8. 
Table 4 (p. 71) gives 2 or 3S values for 2 = 3(2)13 of 6; and 450 such that 
I,(0s) = .95, In(@50) = .50, and also 05/050. 
A. E. 


1197[L].—R. C. T. Smrtu, “Conduction of heat in the semi-infinite solid, 
with a short table of an important integral,’’ Australian Jn. Phys., v. 
6, 1953, p. 127-130. 


Table 1 (p. 128-129) gives 5D values of 


f (1 + u*)— exp [—a(1 + wu’) }du 


for a = .1(.1)2 and U = .1(.1)2, 2.5, 3, ©, and also for a = 2.5, 3, 4, 5 and 
a shorter range of U. 
A. E. 


1198[L].—MicHaEL Trxson, “Tabulation of an integral arising in the 
theory of cooperative phenomena,’’ NBS, Jn. Research, v. 50, 1953, p. 
177-178. 


Table 1. Values of the coefficients c2, in the expansion 
[To(x) P = Lo Com?” 
0 


for m = 0(1)20. Here I(x) is the modified Bessel function of order zero. 
Table 2. Values of 


I(b) = xr f f y [3b — (cos x + cos y + cos z) }'dx dy dz 
0 0 0 


> (2m)! com (36)-2"—2 


to 5D for b-! = .01(.01)1. For 6-' < .8, the expansion in powers of b— was 
used to compute J(6), at most 21 terms of this expansion being required. The 
remaining values of J(b) were obtained by numerical integration. 

A. E. 
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1199[L,V ].—C. TRUESDELL, “Precise theory of the absorption and dispersion 
of forced plane infinitesimal waves according to the Navier-Stokes 
equations,” Jn. Rational Mech. and Analysis, v. 2, 1953, p. 643-741. 
Tables (p. 723-734) computed by Harrison Hancock. Graphs (p. 
735-741). 


These tables give results on the calculation of the propagation of plane 
infinitesimal pressure waves (sound waves) in a uniform fluid governed by 
the Navier-Stokes equations. The fluid is characterized by the viscosity co- 
efficients u and X, the coefficient of heat conduction x, and the specific heats 
Cp, and c,. The tables are divided according to the two parameters, the ratio 
of specific heats 


Y = C/Ce 
and the thermoviscous number Y, 


K 


~ (A + Que," 


Tables given are for (y, Y) of the following values (1, Y) (piezotropic fluids) ; 
(1.10, .5), (1.10, .8), (1.25, .25), (1.25, .5), (1.25, .8), (1.40, .25), (1.40, .5), 


cao), (5,25), (5.8), (5.5), (58), (5105), (5125) 


(1.8, .3), (1.8, .5), (2, .3), (2, .5), (2.2, .3), (2.2, .5) (weak and moderate 
conductors) ; (1.10, 20), (1.10, 30), (1.10, 40), (1.15, 10), (1.15, 20), (1.15, 
25), (1.15, 30), (1.15, 35), (1.15, 40), (1.20, 20), (1.20, 30), (1.20, 40), (1.50, 
50) (strong conductors). In each table, the various dimensionless quantities 
characterizing the speed, the absorption and the dispersion of a plane sound 
wave of circular frequency w are listed against the argument X, the frequency 
number 


Y 


(A + 2) 
Xoo 
YP 
where » the undisturbed pressure of the fluid. The ranges of argument is 
X = .1(.2).7, 1(.5)5. 
H. S. Tsien 


California Institute of Technology 
Pasadena, California 


1200[L ].—G. ZarTARIAN & H. M. Voss, “‘On the evaluation of the function 
fr(M, w),” Jn. Aeron. Sciences, v. 20, 1953, p. 781-782. 


9D table of the coefficients 





n (2M)-4 
Qn(M) = 2 Gram — DI! 
b(M) = 5 — 





joo (F!)*(2(m — 7) + 1)! 
for nm = 0(1)6, M = 5/4, 10/7, 3/2, 5/3, 2, 5/2. These coefficients occur in 
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the expansion 
1 
f e~*#* Jy (wu/M)wdu 
0 
= £(-1) 


n=0 


(- ata) 
N+ 2n+1 “X+2n4+2 , 


and may be used for computing the integral. 
A. E. 


1201[Q].—H. Woon, (a) Kepler’s problem, Roy. Soc. of New South Wales, 
Jn. and Proc., v. 83, 1950, p. 150-163; (b) Kepler’s problem—the 
parabolic case, v. 83, 1950, p. 181-194; (c) Tables for nearly parabolic 
elliptic motion, v. 84, 1951, p. 134-150; (d) Tables for hyperbolic 
motion, v. 84, 1951, p. 151-164; (e) Five figure tables for the calcula- 
tion of ephemerides in parabolic and nearly parabolic motion, Sydney 
Observatory Papers No. 16, 1951. (a), (6), (c) and (d) are also Sydney 
Observatory Papers Nos. 10, 11, 14 and 15 respectively 
The first four papers give the theory and seven-figure tables which are 

especially applicable to the problem of calculating the ephemerides of 

comets in parabolic or nearly parabolic orbits about the sun. The fifth paper 
gives five-figure tables which are sufficiently accurate for finding purposes. 

Some of these tabulations may be of more general mathematical interest. 
The two body problem was solved (kinematically) by Kepler with the 

enunciation of his three laws of planetary motion. The problem of finding the 
coordinates of a planet, or comet, in the plane of its orbit in unperturbed 
motion is Kepler’s problem. If M, EZ and v are the mean, eccentric and true 
anomalies, respectively, and e the eccentricity, then: M = E —e sin E, 
where tan E/2 = ((1 — e)/(1 + ))! tan v/2. This equation is Kepler’s 
equation ; the necessity for its frequent, accurate solution, and the difficulties, 
both numerical and analytical which it presents, have kept alive interest in 
the equation for the past 300 years. In general, e is quite small for the orbits 
of planets, asteroids and most binary stars and there are numerous satis- 
factory methods of solution. Cometary orbits, on the other hand, are very 
often parabolic, or nearly so, and are sometimes even hyperbolic because of 
planetary perturbations, or fictitiously hyperbolic because of observational 
errors. If ¢ is close to unity and E small, special methods must be devised to 
solve Kepler’s equation accurately. Wood writes the equation as follows: 


sine — edu ) 
’ 


D = 12k (1 + e)tg7t = 12¢ + w(1 + 06( ag 


where k is the Gaussian gravitational constant, g the perihelion distance, ¢ 
the time (¢ = 0 at perihelion), « = (1 — e)/(1 + e) and uw = yo/g where xo 
and yo are the rectangular coordinates in the orbital plane with the xo-axis 
directed towards perihelion. 

For e = 0 (e = 1) we have the parabolic case and 


D = 12u = p®, where wp = 2 tan v/2. 


Table 1 in () gives x to 7D for D = 0(0.1)100; to 6D for D = 100(1)1000. 
For D > 1000 use » = Dt — 4/D? + R. Table 2 in (5) gives R to 6D for 
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4/D* = 0(.01)0.41. Previous tables are useful only up to D = 88 and Wood 
states that there have been 36 comets observed outside this range of D, with 
the probability that the future will see an increasingly higher percentage of 
such observations. 

The above tables are also useful in the nearly parabolic case. For e ¥ 1, 
D is redefined as: 

D = 12k(1 + e)ig-ict = 12ce + co', 

and 


w= {7 -AK +R}, 


where c’ and h are certain power series in ¢, and J and K are certain power 
series in eo”. The coefficients of the last terms used in the K and h series have 
been adjusted so that R is negligible in the seventh decimal place. ¢ is 
defined by the top equation and co can be evaluated from Table 1 in (8). 
¢ is tabulated to 7D and & to 5D for e = 0(.001)0.100 in Table 2 in (c) 
(ellipse) and for a(=— e) = 0(.001)0.100 in Table 2 in (d) (hyperbola). 
J and K are given to 7D for e4¢ = 0(.001)0.600 in Table 3 in (c) and sim- 
ilarly for atc in Table 3 in (d). 
Other useful quantities tabulated by Wood are of the form: 


pe sae, I=6 [ice - ner {2 | 
w* wvi—w w* 
where w = e4y for the ellipse and: 
ano{ 2st, r=o | +21, pe nao {Stet}, 
w* wvi+w w* 


where w = ay for the hyperbola. If we define D, as: 


D, = (1 + e) 64 + Ap’, the i = 6(1 + e) + wl. 
A and J are useful in the iterative computation of u and of velocities. The 
other rectangular coordinate in the orbital plane is \ = x0/g = 1 — w®N/2 
(1 + e). This becomes \ = 1 — yu?/4 for the parabolic case. Appropriate 
formulae are given to calculate the equatorial heliocentric coordinates x, y 
and z, from N and uz. 

A and N are tabulated to 7D and J to 4D for e4u = 0(.001)0.600 in 
Table 1 of (c), and similarly for a4u = 0(.001)0.600 in Table 1 of (d). 

Paper (e) gives 5D tables applicable to the three cases considered. There 
has been some slight rearrangement of the formulae and tabulated quantities 
but the details are probably not of sufficient general interest and will not be 
given here. The author states that, in general, the tables have been calcu- 
lated to two extra decimal places. A very useful bibliography will be found 
in paper (qa). 

Joun B. IRWIN 
Goethe Link Observatory 
Indiana University 
Bloomington, Ind. 
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MATHEMATICAL TABLES—ERRATA 


In this issue references have been made to errata in RMT 1165 and 
UMT 184. 


236.—Tosio Kitracawa, Tables of Poisson Distributions. Baifukan, Tokyo, 
1952. 


The above tables have been compared with an unpublished table for 
m = 1(1)10. One misprint and four rounding errors appear as follows: 


m x for read 
7.00 4 .0902262 .0912262 
9.00 15 .0194306 .0194307 
9.00 24 .0000158 .0000159 
10.00 11 -1137363 .1137364 
10.00 15 .0347180 .0347181 


D. TEICHROEW 
NBSINA 


237.—M. Kraitcuik, Recherches sur la Théorie des Nombres, v. 1, Paris, 
1924. 


This work contains (p. 131-191) a table of residue-indices of 2 modulo P 
for P < 300000 as described in UMT 184. 
The following errata are new [for other errata see MTE 107, MTAC, 
v. 2, p. 313]. 
P 


101681 1 
102329 
103183 
103669 
104239 
104383 
104527 
106273 
106663 
106963 
107581 
107941 
108571 
109433 


zg 
8 


ma GO ee te et et 00 © 
CWWFMWNHNN NN WD FS DY 


D. H. L. 


238.—NATIONAL BUREAU OF STANDARDS. Table of Bessel Functions of Frac- 
tional Order, v. 1, New York, Columbia University Press, 1948. 
Function J_;(x) 


x Column For Read 
20.84 & 12305 12299 


WALTER SODEN 
Naval Air Missile Test Center 
Point Mugu, Calif. 
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UNPUBLISHED MATHEMATICAL TABLES 


181[E].—NATIONAL BurEAU oF STANDARDS, Tables of sinh x and cosh x. 
Available on IBM cards at NBSCL. 


A table of sinh x, cosh x for x = 2(.001)10 to 9S has been prepared on 
SEAC. This is an extension of the material available in National Bureau of 
Standards Applied Math. Ser. No. 36. Among those who contributed to the 
preparation of this table were W. F. CAHILL and S. B. Pruscu. 

a Es 


182[E].—NATIONAL BuREAU OF STANDARDS, Table of exp(—x). Available 
on IBM cards at NBSCL. 


A table of exp(—x), x = 2.5(.001)10, to 20D has been prepared on 
SEAC. This is an extension of the material in NBS Applied Math. Ser. No. 
14 which gives exp(—x) for x = — 2.5(.0001)2.5(.001)5(.01)10 to at least 
12D. 

Among those who contributed to the preparation of this table were R. M. 
Davis, E. C. MARDEN, M. PAULSEN and S. B. Pruscu. 

| - 


183[F ].—A. GLopEN, Factorization of 2N* + 1, N < 1000. 17 typewritten 
pages deposited in the UMT FIxe. 


Factorizations are complete up to VN = 100 except for N = 73. They are 
mostly incomplete beyond N = 500. 
A. GLODEN 
rue Jean Jaurés 11 
Luxembourg 


184[F ].—F. GRuENBERGER, Table of residue-indices of 2 for prime moduli 
between 100000 and 110000. One sheet tabulated from punched cards 
deposited in the UMT Fie. Other copies available from the author. 


This table is the beginning of a recalculation of an unreliable table of 
KRaITCHIK! for p < 300000. It gives for each p between the limits indicated 
the integer v = (p — 1)/x where x is the “‘exponent”’ of 2, that is, the least 
positive integer » for which 2” — 1 is divisible by ». New errors uncovered 
in Kraitchik’s table by collation with the present table are given in MTE 
237 in this issue. 

F. GRUENBERGER 
Numerical Analysis Laboratory 
University of Wisconsin 
Madison 6, Wisconsin 


1M. Kraitcaik, Recherches sur la Théorie des Nombres, v. 1, Paris, 1924, p. 131-191. 


185[F].—R. J. Porter, Additional Tables of Irregular Negative Determi- 
nants of Exponent 3n. Typewritten manuscript on deposit in the UMT 
FILE. 

This is an extension of UMT 155 [MTAC, v. 7, p. 34] to determinants 

—D with 50000 < D < 80000. 
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In this range there are three D’s with exponent of irregularity equal to 6, 
namely 
D = 55555, 67899, 70244. 


As in the earlier table, there are many more D’s with exponent 9 in fact the 
following 11. 


52731, 54675, 56403, 58563, 60075, 64395, 70227, 70956, 75411, 76059, 
77571. 
All 1169 other D’s have exponent 3. 

R. J. PORTER 
266 Pickering Road 
Hull, England 


186. L }.—C.-E. Fr6BERG & P. RaBitnow!1tz, Tables of Coulomb Wave Func- 
tions NBS Report, 12 January 1954. Deposited in the UMT Fixe. 

The first volume of a table of Coulomb Wave Functions was issued in the 
NBS Applied Math. Ser. in 1952 (see RMT 1091); it is devoted to the 
regular solution and enables values to about 7 figures to be obtained in the 
range L = 0(1)21, —6< 9 < 6,0< p< 5. Asecond volume will be issued in 
the same series; it will be concerned with the regular and irregular solution 
for L = 0 and in the range 0< n< 10,0< p < 10. This report, of which a 
limited number of copies are available, is concerned witha skeleton table for 
both the regular and irregular solutions. From this table it is possible to 
obtain both the functions to about 6-7 figures in the range L = 0(1)21, 
1 <7 < 10,1 < p < 10, using standard techniques of interpolation and the 
recurrence relations for the functions. 

Although this table itself will be of considerable use to physicists who are 
prepared to do a fairly heavy interpolation job, it is best to regard it as 
the result of proving in a SEAC code for the whole range L = 0(1)25, 
0 <7 < 100,43 < p < 100. From this code, by insertion of suitable param- 
eters, it is possible to obtain a detailed table of the functions in any part of 
the range mentioned. This is prepared, in the first place on magnetic wire, 
and from this IBM cards can be punched, or a Flexowriter manuscript 
obtained. At present it is not thought likely that more detailed tables will 
be printed, but consideration can be given to the preparation of decks or 
manuscripts to satisfy individual or organizational needs. 

A preliminary report on the method of computation, including some new 
integral representations of the functions, was given as Paper No. 9, at the 
NBS Tables Conference on May 15, 1952. 


5. ¥. 


187[L ].— NATIONAL BUREAU OF STANDARDS, Table of the Sievert Integral. 


The Sievert integral is defined by 
6 
y(0,x) = f exp(—-x sec /)dt 


(see Q.19, MTAC, v. 2, p. 196). For 6 = 32, y(@, x) coincides with A7,;(x) 
which has been tabulated by BicKLEY & NAYLER'. 
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This function has been tabulated using SEAC, the computations being 
planned by H. E. SaLzer. Among those who contributed to the preparation 
were R. B. Jasper, P. J. O'HARA, M. PAULSEN, and W. R. SoDERQUIST. 
The quadratures were checked by comparison with the Bickley-Nayler 
table. There is now available, on IBM cards, values to 9S of y(6, x) for 


8 = 0(1°)90°, x = 0(.01)2(.02)5(.05)10. 


A table to about 5S for x = 0(.5)10 and the same values of @ by SIDNEY 
JOHNSTON is referred to in UMT 103, MTAC v. 4, p. 163. 
| A 


1W. G. BickLEy & J. NAYLER, Phil. Mag., s. 7, v. 20, 1935, p. 343-347. 


AUTOMATIC COMPUTING MACHINERY 


Edited by the Staff of the Machine Development Laboratory of the National Bureau of 
Standards. Correspondence regarding the Section should be directed to Dr. E. W. CANNon, 
415 South Building, National Bureau of Standards, Washington 25, D. C. 


TECHNICAL DEVELOPMENTS 
THE DEVELOPMENT OF A.P.E.(X).C. 


1. Historical.—It is the purpose of this paper to give a technical descrip- 
tion of the all purpose electronic digital computer developed at the Birkbeck 
College computation laboratory. Before embarking on this, however, it 
seems appropriate to mention certain historical details which put the present 
machine into perspective. 

In 1946 the present author was engaged in the development of a special 
purpose relay digital calculator for the automatic evaluation of three dimen- 
sional Fourier sums! of the type: 


+H +K +L . —" 
p(x, y,2) = © > |Fiell onl an( = a = + =) — Ahk 
—H —K —-L | a b c 


where the F and a are given in digital form and are available from a punched 
tape. 

The complete design for this machine involved a magnetic disc store for 
256 words of 16 bits, and also a matrix function table (one-many, many-one) 
for producing the cosine function. 

Although this machine was to be relatively compact (600 relays and 100 
vacuum tubes) it did not satisfy aesthetically since the input had to be in 
binary form and output resulted in a like manner. 

Fortunately, at this point, the author was invited to address the Ameri- 
can Society for X-ray and Electron Diffraction, and in the course of this 
visit to the U.S.A. was able to make contact with the pioneer work which was 
being carried out at several centres. 

Through the generosity of the Rockefeller Foundation a period of study 
at the Institute for Advanced Study followed, during which plans were 
completed for an all purpose relay computer using roughly the same equip- 
ment as was already available for the earlier project. 
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Thus, by the end of 1947, the A.R.C. (Automatic Relay Computer) was 
completed. This machine, which is still available for use in this laboratory, 
employs 800 high-speed relays to form a parallel operation single address 
code computer.” The storage organ, which was originally on a nickel plated 
brass drum,’ had a capacity of 256 words each of 21 bit length. This was 
later (vide infra) replaced by an electro-mechanical store for 50 words of the 
same length and the control modified to a pluggable unit for 600 instructions. 

The reason for this seeming volta face lay in the use to which the machine 
was eventually put, which was essentially the Fourier series evaluation 
mentioned earlier. It has proved far more efficient to have a pluggable se- 
quence machine when this is used only intermittently (not more than three 
or four times during a year) and then only by research students. 

Input to the A.R.C. is via teletype tape or push buttons and the output 
is by means of a teleprinter. The code? is as follows: 


No. Symbol Description 
0 TtoM Fill memory from input tape. 
1 TtoTi Start tape moving to 7i and proceed with next order. 
2 Tito M Read material between Ti and 77 into first (7 — 4) positions 
in memory. 
(0) 3 MiytoT Punch contents of memory position 7 to 7 into tape. 
(0) 4 Cto M(x) Shift control to order located at M(x). 
5 Ccto M(x) If number in accumulator > 0 shift control as in 4. 
6 I(n) Shift contents of accumulator and register (A and R) n places 
to left. A(0) A(1) — A(20); R(O) R(1) — R(20) becomes 
A (0), A(2) — A(20), R(O); R(1) — R(20), A (A), etc. 
7 r(n) Shift contents of A n places to right. A(0), A(1) — A(20) 
becomes A (0), A (0), A(1) — A(19), etc. 
8 +MtocA Clear A and add contents of M(x) into it. 
9 +\|M|tocA 
10 -—MtocA 
11 —|M|tocA 
12 +MtoA Add contents of M(x) into present contents of A. 
13 +\|M|toA 
14 -—MtoA 
15 —|M|toA 
16 MRtocA Clear A. Multiply M(x) by R, place 1st 20 digits and sign in 
A after adding unity to 21st. Leave last 20 digits in R. 
17 MRtocA(n.r.o.) As in 16 but without round off. 
18 A+ MtocR Divide A by M(x), place quotient in R with last digit unity. 
Leave remainder in A. 
19 MtocR 
20 RtocA 
21 RtoM 
22 AtoM 
(0) 23 AzrtoM Substitute digits 1-8 of A in order located at M(x). 
24 AtocR 
25 £E Stop and emit end of calculation signal. 





Notice that those instructions marked (0) are no longer available since 
the alteration in control structure, and that the original code-word was 
constructed as follows: 


16 to 20 
Order Code 


9 to 15 
Sequence 


Digit 1 to 8 
Memory 
location (x) 


The zeroth digit was always zero in an instruction and the position 
marked ‘‘sequence’’ was used in orders (2) and (3) to specify (j — 7). 
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The speeds of the various arithmetic operations are: 


+,-. 20 m. sec. 
> ES 1 sec. 
l(n) andr(n) n X 20m. sec. 


In its original form (i.e. on the magnetic drum) transfers to and from M re- 
quired 20 m. sec. (max) ; in the current electro-mechanical device a maximum 
of 250 m. sec. is needed, but this can be reduced by optimum programming, 
to 20 m. sec. 

A partial reason for the removal of the magnetic drum store from A.R.C. 
was a desire to make this unit a part of a small all-electronic calculator. 
This machine, later called S.E.C.* (simple electronic calculator) was never 
brought to a state of useful operation, but served as a test bed for the de- 
signs of the A.P.E.(X).C. machines which followed. The interesting point 
about S.E.C. lies in the small number of vacuum tubes used in its construc- 
tion (230). The code was rudimentary and the input/output took the form of 
switches and neon lamps. A major change in logic was the adoption of a two 
address code® of the type: ‘‘Perform operation A on number in position (x) 
and go to position (y) for the next instruction.” 

This extension removes the chief speed limitation which is inherent in a 
magnetic drum store and in favourable cases, increases the speed of operation 
by a factor of 16-32. 

2. The A.P.E.(X).C. [All Purpose Electronic (X-initial of sponsoring 
agency) Computer |.—At a relatively early stage in the actual construction 
of S.E.C., it was felt that sufficient experience had been gained to justify the 
development of a full-scale machine. A set of standard circuits had been 
devised®:7:* and the problems of drum storage were considered solved. 
Essential details are as follows: 


Drum: Speed 3,800 R.P.M. 
Digit packing 50 bits/inch. 
Track packing 32 tracks/inch. 
Capacity 512 words of 32 bits. 
Repetition rate 30 K.C./sec. 


Tubes: Three major types 6.J.6. 
6.A.L.5. 
6.A.G.5. (only 6 of these). 


Standard Circutts.® 
Binary cell. 
Binary counting cell. 
Two input gate. 
Diode buffer. 
Amplifier inverter. 


No germanium diodes were used, as tests had not shown them to have the 
reliability which was considered needful (in 1949). 

The first machine of A.P.E.(X).C. type was constructed for the British 
Rayon Research Association ; no attempt was made to eliminate redundan- 
cies in the design but the resulting machine, shown in Fig. 1, still had less 
than 500 tubes (no germanium or other semi-conductor diodes). In many 
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Fic. 1. The A.P.E.(R).C. 


ways this machine was still a prototype, the input-output system was not 
fully decided and certain problems connected with track selection had still 
to be solved. Nevertheless from its completion in June 1952, with only 64 
storage positions and teletype input/output, the machine proved reliable 
and useful. 

By January 1953 the complete storage for 512 words was available and 
punched card input/output had been incorporated. From that time, until 
its delivery to the sponsoring agency in July 1953 the machine was in con- 
tinuous use in this laboratory, and it is interesting to note that in the first 
500 hours of operational use only 30 hours of faulty operation were noted 
(94% good operation). Moreover, no regular maintenance was done during 
this period: 


Amongst the problems solved on the above machine (A.P.E.(R).C.) are: 


Octal-Decimal Conversion table 1000 values and differences. 
Decimal-Octal Conversion table 1000 values and differences. 
Tables of binary sines and cosines to 32b. (For machine use). 
Chebyshev polynomial co-efficient evaluation. 

Solution of sets of simultaneous equations (up to 16 X 16). 
Determination of a crystal structure. 
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Perhaps the most interesting of these applications is the last, the structure 
under investigation was that of Oxalic Acid Dihydrate using full three di- 
mensional data. The whole analysis required less than 2 weeks, during which 
time the machine evaluated the equivalent of six sets of three dimensional 
Fourier syntheses (each at 48 points) together with three evaluations of the 
expression :! 


8 ’ o 
F (hkl) = 2 Df, cos an ( + a + a ) 
v=1 a 0 Cc 


each for 1200 values of (hR/). It is estimated that this work would have taken 
a research student between 6 and 12 months by ordinary methods. 

After its move from London to Manchester the machine was in operation 
within one week and it has lately been in use in statistical problems con- 
cerned with mill productivity. 

Whilst A.P.E.(R).C. was being tested several other machines of the 
same general design were constructed. The first group of these were based 
upon a redundancy analysis of A.P.E.(R).C. which resulted in the elimina- 
tion of 180 vacuum tubes from the original system, they are: 


A.P.E.(X).C. at Birkbeck College Computation Laboratory. 

A.P.E.(N).C. at the Norwegian Board for Computing Machines, Oslo. 

M.A.G.I.C. (Magnetic and Germanium Integer Calculator) at 
Wharf Engineering Laboratories, England. 


In the first and last of these machines the storage capacity has been in- 
creased to over 8000 words of 32 bits, and the repetition frequency doubled 
(60 K.C./sec.). Input/output are either by teletype or by punched card. 





Fic. 2. Hollerith Electronic Computer. 2. 
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The second group of machines are based upon the original A.P.E.(R).C. 
design and have been constructed by the British Tabulating Machine Com- 
pany Ltd., they are noticeable for the provision of a special input conversion 
device,’ which enables them to work efficiently in conjunction with decimally 
punched cards. H.E.C.1. (Hollerith Electronic Computer 1) has been in use 
for program investigation for over 12 months and H.E.C.2. was exhibited at 
the Business Efficiency Exhibition, London, in June 1953, giving fault-free 
operation for 992% of the four days of the show. 

A photograph of H.E.C.2., which is now available commercially, is 
shown in Fig. 2. 

3. Code and Control Word Layout on A.P.E.(X).C.—Although minor 
differences both in the actual control word constitution, and in the orders 
themselves, exist between the various members of the A.P.E.(X).C. and 
H.E.C. families ; the scope of these machines will be evident from the follow- 
ing description which refers to A.P.E.(X).C. and M.A.G.I.C. 


11 to 16to 21 to 26 to 
Digits Oto5 6 to 10 15 20 24 25 31 32 
X Track X Loca- Y Track YLoca- Code Spare C6 Spare 
tion tion 


Make up of Instruction Word in A.P.E.(X).C. 


A.P.E.(X).C. Code 


No. Symbol Description 
0 (VE Stop. 
1 (yH Card feed. 
2 (y)P Print first 5 binary digits of number held in register. 
3. (x)(y)Cc If number in A > 0 select next order from position (y) if A < 0 select 
instruction from (x). 
4 (y)l(n) Left shift n places. A and R are connected cyclically. 
5 (y)r(n) Right shift ” places. A (0).A (1)- - -A (31); R(O), R(1)---R(31) becomes 


A (0).A (O)A (1)- - -A (30); A (31), R(O), R(1)- --R(30), ete. 
6 Spare. 


7 Spare. 

8 (x)(y)+c Add contents of (x) into cleared accumulator. 

9 (x)(y)-—c Subtract contents of (x) from cleared accumulator. 
10 (x)(y)+ Add contents of (x) into accumulator. 
11 (x)(y)- Subtract contents of (x) from accumulator. 


12 (x)(y)+R Transfer contents of (x) to R. 

13. (x)(v)X(m) Multiply number in (x) by last » digits of number in R. First half of 
product in A last part in R. 

14. (x)(y) A(z) Record first (or last) » digits of A in (x). 

15 (x)(y) R(m) Record first (or last) m digits of R in (x). 


Some features of this code require further description. In the first place 
it will be noticed that each order has both an (x) and a (y) address. The first 
(x) refers to the location of the number required in the performance of the 
operation, whilst the second (y) gives the location of the next instruction. 
The number () which specifies shift, precision of multiplication, and num- 
bers of digits to be recorded has to be placed in digits 26-31 (Marked C6). 
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It is chosen in accord with the following rules: 


l(n) C6 contains (7) 

r(n) C6 contains (64-7) 

X (n) C6 contains (64-7) 

A(n) C6 contains (32-m) if 1st (least significant) 1 
digits of A are to be re- 
corded, but (64-2) if 
first n digits are to be 
recorded. 

R(n) C6 as in A(n) 


The multiplier is based upon the non-restoring method described elsewhere,'” 
and operates upon numbers of any sign; it normally short-cuts and the pro- 
vision of the () digit facility is an aid to programming rather than to speed. 

The times of the various arithmetic operations (in A.P.E.(X).C.) are: 


+,—, 600 sec. 
xX m X 600u sec. 
r(n) 600u sec. 
I(n) 1,200u sec. 


In these figures no account has been taken of access; with efficient optimum 
programming it is found that an increase by a factor of only slightly greater 
than two is appropriate, but it must be remembered that by “‘non-optimum”’ 
programming all of these figures can be increased by a factor of up to 32! 
It will be seen that the multiplication time is given as m X 600y sec. The 
factor m represents the number of changes 0 to 1 or 1 to 0 of the digits in the 
multiplier. 

A point worthy of mention is that the correct working of the machine 
in no way depends upon the accuracy of the optimum programming, if the 
coder makes a mistake in this it will not cause a machine fault but merely a 
waste of time. This is a feature which is absent in at least one other machine 
employing a multi-address code. To complete the description, it must be 
mentioned that each track on the drum contains 32 words, and it is within 
these that optimisation is required since inter-track switching is performed 
between words and causes no delay. 

The storage (8 X 32 tracks each of 32 words) is at present available only 
in blocks of 32 tracks and the programmer has to select and plug those 
tracks which he will require in a given problem. This limitation may or may 
not be temporary, it is in contemplation to have a telephone director type 
system which will perform the track selection upon instruction either from 
the machine, or from punched cards; the building of this is, however, await- 
ing the result of coding studies which are still in progress. 

ANDREW D. BootH 
Birkbeck College Computation Lab. 
University of London 
England 


1A. D. Booru, Fourier re in X-ray Organic Structure Analysis. Cambridge, 1948. 

2A. D. Boor & K. H. V. BritTE N, Coding for A.R.C. . Mimeographed. Princeton, a 

3A. D. Boorn, “A betb 4 digital storage system,” Electronic Eng., v. 21, 1949, 
234-238. 
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4N. Kitz, Thesis for the degree of M.Sc. University of London, 1951. 

5A. D. Boota & K. H. V. Bootu, Automatic Digital Calculators, London, 1953. 

6 A. D. Boor, “The physical realization of an electronic digital computer,” Electronic 
Eng., v. 22, 1950, p. 492-498. 

7A. D. Booru, “The physical realization of electronic digital computer,” Part II. 
Electronic Eng., v. 24, 1952, p. 442-445. 

8 A. D. Boots, “On optimum relations between circuit elements and logical symbols in 
design of electronic calculators.” Brit. I.R.E. Jn., v. 12, 1952, p. 587-594. 
py ° Birp, “Computing machines—Input and output,” Electronic Eng., v. 25, 1953, p. 

#” A. D. Booru, “A signed binary multiplication technique,” Q. Jn. Mech. Appl. Math., 
v. 4, 1941, p. 236-240. 


BIBLIOGRAPHY OF CODING PROCEDURE 


The material described below is among that which has been added to the 
collection at the National Bureau of Standards Computation Laboratory. 
A similar collection is available at the National Bureau of Standards’ 
Institute for Numerical Analysis. 

Material for inclusion in these collections should be sent to the Computa- 
tion Laboratory, National Bureau of Standards, Washington 25, D. C., for 
the attention of J. H. WEGSTEIN. 


35. UNIVERSITY OF ToRONTO, FERUT, Program Notes, Library Routines. 
This is an expandable library of FERUT routines preceded by 
chapters from the Programmers’ Handbook for the Manchester Electronic 

Computer and notes on the FERUT library. 

FERRANTI Ltp., Introduction to programming, notes, and routine speci- 

fications for the Manchester Electronic Digital Computer. 

This booklet was prepared as an introduction to programming to 
facilitate the study of the Programmers’ Handbook for the Manchester 
Computer Mark IT. It is prepared in loose-leaf book form so that Ferranti 
Routine Specifications can be added as they are issued. 

37. THE CompuTING MACHINE LABORATORY, THE UNIVERSITY OF MAN- 
CHESTER, Programmers’ Handbook (2nd Edition) for the Manchester 
Electronic Computer Mark II. 

This contains chapters on the logical design of the machine and the 
instruction code, coding examples, the preparation of a problem with 
many programming details, aids to coding (e.g., translation routines, 
interpretive routines), the calculation of functions of a single variable, 
numerical solution of ordinary differential equations, fault diagnosis in 
programmes, measures against machine breakdown, and the console. 

38. KLAUS JAMELSON & FREDERICH L. BAvUER, “‘Optimale Rechengenauig- 
keit bei Rechenanlagen mit gleitendem Komma,”’ Z. angew. Math. Phys. 
v. 4, 1953, p. 312-316. 

In most machines with floating decimal (or binary) point a complete 
“zero-shift’’ or “normalization” is carried out either immediately after 
each operation (Mark II, Bell-computer) or before a multiplication 
(BARK), i.e., the numbers are shifted to the left until the first non- 
vanishing digit occupies the leftmost position. This method has the 
drawback that occasionally (e.g., after a subtraction of two nearly equal 
numbers) insignificant figures, due to round-off etc., are shifted far to 
the left, where they may spoil subsequent operations. The paper under 
review proposes to make incomplete shifts in such cases, keeping the 


36 
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insignificant figures at the right. The allowable amount of shift and how 
it is determined by the machine itself is discussed. For binary computers, 
since the first digit of a mantissa is now not necessarily different from 
zero, a modification of the usual conversion routines becomes necessary, 
and this is discussed in detail. 

PETER HENRICI 


NBSCL 


39. 


40. 


41. 


T. Pearcy (Commonwealth Scientific and Industrial Research Organ- 
ization (CSIRO)), Automatic Computation Part I, The Design of the 
Mk. I Automatic Computer. 

This is a description including logical block diagrams of the design of 
the Mk. I computer. it is a 1024 word, mercury delay-line memory 
machine. 

T. Pearcey (C.S.1.R.O.), Automatic Computation Part II, Programmes 
for an Automatic Computer. 

This is an introductory programmer’s manual for the CSIRO Mk. I 
computer. 

T. Pearcey (C.S.1.R.O.), Automatic Computation Part III, Programmes 
for the Mk. I Computer: Pt. I. 

This is an extension of Part II and discusses the preparation of sub- 
routines for the Mk. I computer. 

IBM ScIentiFIC CoMPUTING SERVICE, IBM 701 Speedcoding System. 

The ‘“‘speedcoding system”’ refers to the IBM 701 equipment and a 
method in which the computer aids in the preparation of its own pro- 
grams. The system permits the programmer to supply the binary ma- 
chine with input data in decimal form and with a pseudo-code. The 
pseudo-code is more compact and can be more easily and accurately 
written than the true machine code. The objective of this method is to 
shorten the time required in preparing problems for solution with the 
high speed computer. 

Speedcoding utilizes a floating point system. The pseudo-code in- 
cludes the regular instructions of the true code but also includes single 
instructions which designate computation of square root, evaluation of 
trigonometric and exponential functions, certain tape and drum manipu- 
lations, and error checking operations. 

This booklet is a manual for users of the speedcoding system. 

J. H. WEGSsTEIN 


NBSCL 


43. 


R. R. Hare, Coding for the Florida Automatic Computer (FLAC) (re- 
vised April 1953). 

This report is intended as a reference for the coding of problems for 
the Florida Automatic Computer at Patrick Air Force Base, Florida. 
It briefly describes the floating three-address system employed by the 
FLAC, together with the instructions which may be performed by this 
computer. 


. MATHEMATICS PANEL, OAK RIDGE NATIONAL LABORATORY, Manual for 


the ORACLE. 





45. 
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This beautifully printed and handsomely bound manual is expand- 

able. It presently contains the following sections: “Electronic Digital 
Computers” by C. L. PERRY, a general description of digital computers; 
“‘Basic Operations of the ORACLE” by A. S. HousEHo.per, which in- 
cludes a list of the ORACLE operations; ‘Flow Charts and Coding”’ by 
J. Mosman, which explains flow diagramming and gives examples of 
flow diagrams and their corresponding codes. 
H. RUTISHAUSER, Automatische Rechenplanfertigung bei programm- 
gesteuerten Rechenmaschinen. Mitteilungen aus dem Institut fiir ange- 
wandte Mathematik an der Eidgenéssischen Technischen Hochschule 
in Ziirich, No. 3. 45 p., Basel 1952. 

A considerable part of the time consumed by the solution of a prob- 
lem on an automatic computing machine is usually spent for coding. 
Various attempts have therefore been made to mechanise not only the 
computation, but also the coding procedure (Coding machine for Mark 
III, K. Zuse’s “Plankalkiil” etc.). All these methods require special 
equipment, and the computing machine is not used for the coding. Con- 
trary to this, in the report under review H. RUTISHAUSER gives a 
thorough and complete description of his method to “compute” the code 
for a given problem on the computing machine itself. Although this 
method is in principle quite general and could, with some modifications, 
be carried out on any automatic computing machine, the present exposi- 
tion is based upon the special properties of the machine to be built at the 
Institute for Applied Mathematics in Ziirich. 

The first section of the report is devoted to a short description of this 
machine, the plans of which have grown out of the experience of its 
planners with Zuse’s and Aiken’s machines. It is a one-address machine 
with an internal drum memory of 1000 cells. Instructions are stored in a 
separate memory of the same capacity, each cell holding two instruc- 
tions. Orders can be modified by means of nine index registers. A special 
order is provided for the punching of an order computed on the machine. 

The second section is the heart of the paper. It describes the computa- 
tion of a code for the evaluation of a given algebraic expression (“bracket 
expression’’). Associated with every symbol EF; of a bracket expression 
(brackets and operation symbols included) are two integers, a, and ,. 
Roughly speaking, the more in the interior a symbol is located in the 
bracket expression, the larger is the corresponding number a;. By an 
ingenious algorithm the bracket expression is then decomposed step by 
step, starting from the interior-most elements, and the corresponding 
orders are computed. 

The third section deals with cyclic problems. Here one and the same 
bracket expression is to be evaluated several times and for several values 
(and/or addresses) of the numerical parameters. The author proposes to 
“‘stretch”’ the code, i.e., to unfold the cycles and to compute and to store 
the instructions for each iteration step separately. It is true that some 
computation time is thus saved by avoiding a decision as to the correct 
alternative after completion of each cycle. In view of the obvious limita- 
tions of memory space, the procedure nevertheless does not seem to be 
quite convincing. 
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In the fourth and last section an example is given for the computation 
of a code for a cyclic problem without unfolding. Here the general method 
is complicated by a considerable amount of special rules and provisions. 
Only experience can show if the author’s method results in a real gain 
over straight-forward coding by hand. 

In the whole this is a remarkable publication, which should provide 
a sound basis for further researches in the theory of automatic coding. 

PETER HENRICI 


NBSCL 


46. 


47. 


49. 


National Bureau of Standards Report, September 30, 1952. SEAC 

Operating and Programming Notes, V. 

26. Interpretive Subroutine for Operations on Complex Numbers 

27. Subroutine for Square Root of a Single Precision Number with 
Floating Binary Point 

28. Subroutine for Sinh y, Cosh y, |y| < 2.0634 

29. Code Checking 

30. Basic Arithmetic Operations for Double Precision Numbers with 
Fixed Binary Point 

31. Interpretive Subroutine for Operations on Double Precision Num- 
bers 

National Bureau of Standards Report, April 13, 1953. SEAC Operating 

and Programming Notes, V1. 

32. Interpretive Subroutines $36 and S37 for Operations on Single Pre- 
cision Numbers in Floating Binary Point Form, and Breakpoint 
Subroutine S38 

33. Subroutine for Double Precision, Floating Decimal Point Operations 

34. Subroutine for Log.N; Single Precision, Floating Binary Point 

35. Subroutine Assembly Routine (SASS) 

36. Tape Control Subroutines $21 and S22 


. National Bureau of Standards Report, April 24, 1953. SEAC Operating 


and Programming Notes, VII. 

37. Basic Arithmetic Operations II, Floating Binary Point, Single Pre- 
cision Numbers 

38. Interpretive Subroutine for Operations on Double Precision Num- 
bers with Floating Decimal Point 

39. Subroutine for the Square Root of a Double Precision Number with 
Fixed Binary Point 

40. Subroutine for the Square Root of a Single Precision Number with 
Floating Binary Point II 

41. Automonitor Processing Routine (BIG BOIE) 

National Bureau of Standards Report, November 13, 1953. SEAC 

Operating and Programming Notes, VIII. 

42. Subroutines S30, S31, S32, and S33 for the step-by-step integration 
of the first order differential equation y’ = f(x, y) 

43. Subroutine for 16-point Gaussian quadrature 

44. Subroutine for binary to decimal conversion of a double precision 
number with fixed binary point 





50. 


51. 


NB 


52. 
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45. Subroutine for binary to decimal conversion of a single precision 
number with floating binary point 
46. Subroutine for decimal to binary conversion of a single precision 
number with floating binary point 
47. Subroutine for conversion of decimal degrees to binary radians and 
vice versa 
48. Subroutine for the evaluation of polynomials, single precision, fixed 
binary point 
49. Automatic transfer routine 
50. NBS INsTITUTE FOR NUMERICAL ANALYsiIs. SWAC Memoranda. 
SWAC Memoranda are prepared in order to maintain a record of 
current developments in SWAC operating and programming. These 
memoranda include routines of general interest and applicability, facts 
of interest to programmers concerning modifications or newly installed 
facilities of the SWAC, and other information of general importance. 
The information incorporated in these memoranda represent the efforts 
being made in the Institute to develop techniques (a) for reducing the 
amount of time that goes into code preparation and checking, (b) for 
checking machine and auxiliary equipment operations, (c) for mechaniz- 
ing as much as possible the code preparation process, and (d) for making 
optimum use of SWAC equipment. 
Memoranda in this series will in general be reviewed by title only. 
2. Interpretation Routine 
3. General Test Routine 
4. The Magnetic Drum 
5. Break Point Digit 
6. Relative (Symbolic) Coding 
7 
8. 
L 


. Diagnostic Extract Routine 
Memory Test Series 00000 to 00099 
. J. Paice & C. B. Tompkins, Systematic generation of permutations on 
an automatic computer and an application to a problem concerning finite 
groups. NBS working paper, January 30, 1953. 

The authors give a method for generating permutations on m marks 
which they have coded for SWAC for n = 7,9. The method is as follows: 
Suppose that all permutations on m — 1 marks are known. If (41, 42, 
+++, dni) is such a permutation, then (7, i2,---, ins, # — 1) and all its 
cyclic permutations will generate all permutations on m marks. The 
method is applied to the problem of determining all ‘“‘complete’’ mappings 
of a certain group of order 8 onto itself. An alternative method of gener- 
ating permutations due to MARSHALL HALL is discussed. A detailed 
code is included. 


51. 


Morris NEWMAN 
NBSCL 


52. ACM, Proceedings of the meeting at Pittsburgh, Pa., May 2-3, 1952. 
Of the 41 papers in these proceedings the following are noted under 
this bibliography: 
“Construction and Use of Subroutines for the SEAC,” J. H. Levin. 
“The Use of Subroutines on SWAC,” ROSELYN LIPKIS. 











53. 


110 


55. 
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“The Use of Subroutines in Programmes,’”’ Davip J. WHEELER. 
“Progress of the Whirlwind Computer Towards an Automatic Pro- 
gramming Procedure,” JOHN W. Carr. 
“The Education of a Computer,’’ Grace M. Hopper. 
ACM, Proceedings of the meeting at Toronto, Ontario, September 8-10, 
1952. 

Of the 35 papers in these proceedings the following are noted under 
this bibliography: 
“Compiling Routines,” R. K. RipGway. 
‘Machine Aids to Coding,” E. J. ISAac. 
“Computer Aids to Code Checking,” I. C. DrExM. 
“Input Scaling and Output Scaling for a Binary Calculator,” E. F. 
Copp, H. L. HERRICK. 
“Logical or Non-Mathematical Programmes,” C. S. STRACHEY. 
“A Simplified Universal Turing Machine,”’ E. F. Moore. 
“Simple Learning by a Digital Computer,” The CompuTaTIon LABorRA- 
Tory, Harvard University. 
“Use of Continued Fractions in High-Speed Computing,’’ D. TErcu- 
ROEW. 
“Interpretive Sub-Routines,” J. M. Bennett, D. G. Prinz, M. L. 
Woops. 
“The Use of Sub-Routines on SEAC for Numerical Integrations of 
Differential Equations and for Gaussian Quadrature,’’ P. RaABINOWITZ. 
“Pure and Applied Programming,” M. V. WILKEs. 


. NATIONAL PuHysIcAL LABORATORY, Teddington, England. Automatic 


Digital Computation Symposium, March 25-28, 1953. 

This is a collection of copies of most of the 37 papers presented. The 

following are noted under this bibliography: 

“Optimum Coding,” G. G. ALway. 

“Micro Programming and the Choice of Order Code,” J. G. STRINGER. 
“Conversion Routines,’’ E. N. Mutcu. 

“Getting Programmes Right,’’ S. GILL. 

“Diagnostic Programmes,” R. L. GRIMSDALE. 

OsBornE, E. E., Solution of the Matrix Equation (M — QD)X = 0. 
NBS report, October 16, 1953. 

The author gives a solution by the power method (successive 
matrix-vector multiplication and vector normalization) of the equation 
of the title, M being a 6 X 6 matrix having complex elements and D a 
non-singular real diagonal matrix. Setting A = D~'M, the equation be- 
comes the more familiar one 


(A —QNX = 0, 


and this is the one the author actually solves. 

The preliminary multiplication D-'M was performed on IBM equip- 
ment. The computation of the eigenvalues Q and eigenvectors X was 
programmed for SWAC. Detailed flow charts and codes are given for 
two such programs, one which can be stored in the high speed memory, 
the other requiring external storage capacity such as a magnetic drum. 

Morris NEWMAN 
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BIBLIOGRAPHY Z 


1114. Anon., ‘‘New knowledge—computers are busy these days,” Westing- 

house Engineer, Jan. 1954, p. 62-63. 

Industry is finding computers economical as well as efficient. An Eastern 
utility is using one to establish power-system loss formulas and data for load 
dispatching on large interconnected systems at a saving of $60,000 yearly. 

They are being used for the automatic solution of problems involved in 
the routine application of electrical machinery ; problems in connection with 
blades of compressors and steam turbines; reduction of data on all effects of 
inner-cooled generators and stress-analysis problems on rigid steam piping. 

A. R. Cock 
NBSCL 


1115. ANnon., ‘‘New digital computer speeds math solution at Oak Ridge,” 

Industrial Laboratories, Nov. 1953, p. 128-129. 

The Oak Ridge Automatic Computer Logical Engine, ORACLE, is 
described in 9 paragraphs and 3 pictures. It is a descendent of the Institute 
for Advanced Studies Computer, is fully parallel, and has cathode-ray tubes 
storage for 2048 words. Net addition time is 5 microseconds and multiplica- 
tion time 500 microseconds. 

R. D. ELBouRN 
NBSCML 


1116. ANon., “Utilization of digital computing machines,’’ Nature, v. 172, 

1953, p. 649-651. 

This is a report of a meeting of the Engineering Section of the British 
Association in September 1953. The speakers were F. C. Wriitams, R. K. 
LIVESLEY and G. G. ALway. It is interesting to note British experience and 
opinions in this field: for instance, that at present the need is for more pro- 
grammers and programs rather than for more machines, that it takes about 
three months to train a programmer and about the same time to prepare a 
significant problem, and that the speed factor is 150: 1 over hand computing. 
The desirability of making the machine carry out as much as possible of the 
preparation of data, no matter how trivial it is, is stressed. It is also noted 
that the knowledge of the capabilities of computers gained during the solu- 
tion of one problem has often led to customers realization of the feasibility of 
attacking other problems in their field. 

Among the problems which have been handled (on the machines at 
University of Manchester and at the National Physical Laboratory) are: 
torsional and lateral oscillations of systems of shafts, design of cams, 
analysis of rigid-frame structures, flutter analysis. 

> 


1117. P. R. BAGLEy, Electronic Digital Machines for High-speed Information 
Searching. Master’s Dissertation, Mass. Inst. of Tech., Cambridge, 1951, 
ix + 133 p. 
The author’s objective is the examination of the suitability of methods 
and machines for the high-speed location of related items in a large body of 
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indexed information. His study is a timely one because of the mountainous 
volume of information, growing at an ever-increasing rate, which exists in 
the scientific and engineering fields. For example, a search of Chemical 
Abstracts today involves over a million abstracts, with the present trend of 
publication indicating that the total abstracts published in this field by 1960 
will be almost 1,800,000. In another representative field, Electrical Engineer- 
ing, it has been estimated that by 1960 the total of published papers of inter- 
est may reach 700,000. 

Since the effective planning of research requires a knowledge of previous 
activity in the field, it is apparent that the increasing size of our ‘‘storehouse”’ 
of information makes it imperative that we attempt to increase the speed of 
our information searching facility. A notable stride in this direction has been 
taken in the development of the Rapid Selector, designed by Dr. VANNEVAR 
Bus3, which is in use and being developed further, by the U. S. Department 
of Agriculture Library. In respect to the problem of information searching, 
methods and techniques for searching inventory records are interesting. For 
example, in the October 1953, Proceedings of the I.R.E., there was reported 
an investigation of photographic techniques for information storage under- 
way by the International Telemeter Corporation, which could result in 
further alleviation of the information searching problem. 

The thesis is divided into three parts: Part One, Introduction to the 
Information Searching Problem; Part Two, Investigation of Electronic 
Digital Methods for High-Speed Scanning and Selection; Part Three, Con- 
clusions and Recommendations. Four Appendixes are added: Appendix 1, 
Derivation of Superposed Coding Formulas; Appendix 2, Description of 
Major Electronic Computer Components; Appendix 3, Guide to Coding for 
the Whirlwind I Computer; Appendix 4, Searching with the Whirlwind I 
Computer. A thorough analysis of the problem prefaces an attempt at solv- 
ing it by use of a high-speed electronic digital computer. 

The conclusion of the author is that the general-purpose digital computer 
is virtually disqualified as a selection device because of its sequential opera- 
tion and the fact that, in the application, substantial time is required for the 
transfer of numbers back and forth between the storage element and the 
arithmetic element of the machine. He considers it feasible to develop an 
electronic computer more suitable for information searching than the general- 
purpose digital computer, and presents a block diagram for such a device. 

E. W. C. 


1118. H. Buckner, F. J. Weyt, L. Brermann & K. ZuseE, Probleme der 
Entwicklung programmgesteuerter Rechengeraite und Integrieranlagen. 
Rhein.-Westf. Technische Hochschule Aachen, Mathematisches Institut, 
H. CrEMER, ed., Aachen, 1953, xiii + 75 p. 

This book contains the essential contents of four speeches given at a col- 
loquium in Aachen in July 1952, under the sponsorship of the Institute for 
Mathematics, Mechanics, Physics and Theoretical Physics of the Rheinish- 
Westfalischen Technischen Hochschule. 

The four chapters of the book, based upon the speeches of the authors in 
the order above, are the following: Uber die Entwicklung des Integromat; 
Aufbauprinzip, Arbeitsweise und Leistungsfahigkeit elektronischer pro- 
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grammgesteuerter Rechenautomaten und ihre Bedeutung fiir die naturwis- 
senschaftliche Forschung; Die Gottinger Entwicklungen elektronischer 
Rechenautomaten ; Uber programmgesteuerte Rechengerate fiir industrielle 
Verwendung. 

In the first chapter, BUCKNER gives a functional description of the 
Integromat, a differential analyzer. In the second chapter, by F. J. WeyL, 
the characteristics of electronic digital computers now in use in the United 
States are compared and some of their uses are discussed. Included in the 
comparison of machines is the primary motivation of their construction— 
ballistics computations, the development of improved components and more 
flexible computational systems, or commercial and industrial applications. 
In the third chapter, L. BrERMANN discusses the features of a new digital 
computer in operation at the Max-Planck Institute of Physics. The com- 
puter is a moderate-speed, tape-fed, magnetic drum machine, operating on 
32 binary-digit numbers, with a 3 binary-digit integral part. Some of the 
problems being solved on the computer are described. The last chapter, by 
ZUSE, is an exposition of his development of relay computers, prefaced by a 
review of the difficulties he encountered up to the time of the installation of a 
Zuse computer in the Eidgendssische Technische Hochschule at Ziirich in 
1950. He presents an interesting explanation of the reasoning underlying 
many of his design decisions, and states a case for moderate speed, simple- 
to-operate computers, particularly for application to engineering problems. 

The book closes with a recording of the discussion which ensued after the 
presentation of the papers. 


E. W. C. 


1119. OrFIcE oF NAVAL RESEARCH, A Survey of Automatic Digital Com- 
puters. Washington, 1953, vi + 109 p. Available as PB111293, U. S. 
Dept. of Commerce. Office of Technical Services, Washington 25, D. C. 
Price $2.00. 

This compendium of data on automatic digital computers incorporates 
the results of surveys by the Flight Research Laboratory, Wright-Patterson 
Air Force Base, and the Mathematical Sciences Division, Office of Naval 
Research. 

Included in the survey is one page of information on each of ninety-eight 
computers, among which are included machines in use, or under development, 
in Australia, Belgium, Canada, England, France, Germany, Holland, Japan 
Norway, Sweden, Switzerland and the United States. 

A partial listing of the information given on the computers included in 
the survey is the builder, location of installations, availability of program- 
ming service and computing time, operating schedule, number base, word 
length, instruction type, sequence control, built-in operations and time re- 
quired for their execution, description of internal storage and listing of in- 
put-output devices. 

The Office of Naval Research survey is a useful reference for workers in 
the automatic digital computer field and also for those who, though not con- 
cerned with their development, are interested in the application of high- 
speed digital computers to their problems. 

E. W. C. 
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1120. Mina REEs, “Applying Computers to Machine Control,’’ Machine 
Design, v. 25, 1953, p. 324-334. 
An abstract of an address entitled ‘‘The Future Use of Digital Com- 
puters” given at Eighteenth National Applied Mechanics Conference of the 
American Society of Mechanical Engineers, Minneapolis, June 1953. 


NEws 


Eastern Joint Computer Conference and Exhibition —The Joint Computer Conference 
on Information Processing Systems—Reliability and Requirements was held at the Statler 
Hotel in Washington, D. C., December 8-10, 1953, under the joint sponsorship of the 
American Institute of Electrical Engineers, the Institute of Radio Engineers, and the Associ- 





tion for Computing Machinery. 
Session I Tuesday, December 8, 10:15 AM 


Chairman—F. J. MAGInniss, Secretary, AIEE Committee on Computing Devices 


Address of Welcome 
Keynote Address 


The RTMA Support of the 1950 Computer 
Conference—A Progress Report 

The Use of Electronic Data Processing Sys- 
tems in the Life Insurance Business 

Computer Applications in Air Traffic Con- 
trol 

Discussion of papers 

Session II Tuesday, December 8, 2:00 PM 


Joun H. Howarp, Chairman, IRE Pro- 
fessional Group, Electronic Computers 
Howarp T. ENGstroM, Remington Rand, 

Inc. 
Tuomas H. Briccs, Burroughs Corporation 


M. E. Davis, Metropolitan Life Insurance 
Company 

VERNON I. WEIBE, Air Transport Associa- 
tion of America 


Chairman—Howarp T. ENGstrom, Remington Rand, Inc. 


Data Processing Requirements for the Pur- 
poses of Numerical Weather Prediction 
Methods Used to Improve Reliability in 
Military Electronics Equipment 

Digital Computers for Linear, Real-Time 
Control Systems 

The MIT Magnetic-Core Memory 


Discussion of papers 


JosEPH SMAGORINSKY, U. S. WEATHER 
Bureau 
L. D. WuITELock, Bureau of Ships 


Rapu B. Conn, Jet Prop. Lab., California 
Institute of Technology 

WILLIAM T. PAPIAN, Massachusetts Insti- 
tute of Technology 


Session III Wednesday, December 9, 9:00 AM 
Chairman—C. V. L. Smitu, Office of Naval Research 


Reliability Experience on the OARAC 


Operating Experience with the Los Alamos 
701 

Acceptance Tests for the Raytheon Hurri- 
cane Computer 

Reliability of a Large REAC(R) Installa- 
tion 

National Bureau of Standards Performance 
Tests 

Experience on the Air Force UNIVAC 


Discussion of papers 


Session IV Wednesday, December 9, 2:00 PM 
Chairman—C. C. BRaMBLE, U. S. Naval Proving Ground, Dahlgren 


RoBert W. House, Wright-Patterson Air 
Force Base 

Wriu1aM G. Bouricrus, Los Alamos Scien- 
tific Laboratory 

Francis J. Murray, Columbia University 


BERNARD D. LoveMAN, Reeves Instrument 
Corporation 

S. N. ALEXANDER & R. D. ELsourn, 
National Bureau of Standards 

RoBert Kopp, Headquaters, U. S. Air 
Force 
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Electron Tube and Crystal Diode Experi- 
ence in Computing Equipment 

Reliability and Characteristics of 
ILLIAC Electrostatic Memory 

Electron Tube Performance in Some Typical 
Military Environments 

Discussion of papers 


the 


Session V Thursday, December 10, 9:00 AM 


115 
J. A. Goetz & H. J. Geiser, IBM Corpora- 
tion 


JoserH M. Wier, University of Iillinois 


D. W. SHarp, Aeronautical Radio, Inc. 


Chairman—R. F. CLIPPINGER, Raytheon Manufacturing Company 


SEAC—Review of Three Years of Opera- 
tion 

A Review of ORDVAC Operating Experi- 
ence 

Some Remarks on Logical Design and Pro- 
gramming Checks 

The Advantages of Built-in Checking 

Recent Progress in the Production of Error 
Free Magnetic Computer Tape 

Discussion of papers 


Session VI Thursday, December 10, 2:00 PM 
Chairman—S. B. Williams, President, ACM 


Reliability of Electrolytic Capacitors in 
Computers 

Reliability and its Relation to Suitability 
and Predictability 

Case Histories in Resistor Stability 


Discussion of papers 
Summary 


P. D. Saupe, Jr., & R. A. Krrscu, National 
Bureau of Standards 

CuHarLes R. WitiiaMs, Ballistic Research 
Laboratory 

HERMAN H. GotpsTINE, The Institute for 
Advanced Study 

Joun W. Mauca ty, Remington Rand, Inc. 

J. C. Caapman & W. W. WETZEL, Minne- 
sota Mining and Manufacturing Company 


Mark VAN Buskirk, P. R. MALLory and 
Company, Inc. 
E. B. FERRELL, Bell Telephone Laboratories 


JEssE MArsTEN, International Resistance 
Company 


ALLEN V. Ast1n, Director, National Bureau 
of Standards 


Group discussions were held each day at 4 PM which included the following topics: 


Wednesday, December 9 
Applications—Technical 


Applications—Commercial and Industrial 


Diagnostic Checks 
Magnetic Tape Standards 


Crystal Diodes (design for reliability) 


GEORGE PETRIE, Chairman, International 
Business Machines Corp. 

WALTER FRESE, Chairman, General Ac- 
counting Office 

J. J. Eacuus, Chairman, Dept. of Defense 

ALLEN SHouP, Chairman, Shoup Engineer- 
ing Co. 

Rapa J. Stutz, Chairman, National Bureau 
of Standards 
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1121. W. G. ANDERSON & E. H. FRiTzeE, “Instrument approach system 

steering computer,” I. R. E., Proc., v. 41, 1953, p. 219-228. 

This paper is mainly devoted to a theoretical discussion of a computer to 
be used in providing steering indications to an aircraft pilot during a blind 
landing using the Instrument Landing Approach System. The purpose of 
the computer considered is to combine signals from three sources used by a 
pilot in flying such an approach. These sources are the vertical gyro (bank 
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angle), the directional gyro (deviation of heading from runway direction) 
and localizer receiver (lateral deviation from desired flight path). 

The approach procedure is interpreted as a feedback control process. 
The characteristics of the data sources and the aerodynamic characteristics 
of the aircraft are considered in order to devise a steering indication capable 
of yielding suitable approaches. In particular, a system for deriving a desir- 
able signal for rate of change of lateral deviation is described. 

The computer is described briefly. It is of conventional 400-cycle analog 
type. The use of vacuum tubes is held to a minimum and limiters are em- 
ployed to prevent indication of excessive bank angles. 

L. H. O'NEILL 
Columbia University 
New York, New York 


1122. A. R. BootHroyp, “Design of electric wave filters with the aid of the 
electrolytic tank,” Inst. Elect. Eng., Proc., v. 98, part IV, 1951, p. 65-93. 
This article is intended as a ‘‘treatise”’ on the subject. An appendix gives 

details of tank construction. 


1123. JoHN BROOMALL & LEON RIEBMAN, “‘A sampling analogue computer,” 

I. R. E., Proc., v. 40, 1952, p. 568-572. 

The analogue computer described in this article is based on an idea of 
C. J. Herscu and J. F. FELKer. The input of the computer consists of three 
steady d. c. voltages with values X, Y and Z. The computation is based on 
the sampling of d. c. voltages by means of electronic switching circuits using 
diodes. The computer has two units; the first of these is an algebraic unit, 
which produces a voltage pulse whose amplitude W is given by W = RZ YX—. 
The second unit produces a steady voltage U in the form aZ YX— + 6b for 
constants a and b. In the algebraic unit the voltages Z and X are sampled 
and the condensers of identical RC networks are charged to these voltages. 
At a time ¢ after this sampling, the two condensers have voltages Xexp 
(— t/RC) and Zexp(— t/RC). The first of these voltages is compared by an 
electronic circuit with the voltage Y and when these are equal a sample pulse 
whose size is proportional to the second voltage is emitted by this unit. At 
the time to when the first of these voltages equals Y, exp(— to/RC) = Y/X, 
and hence at this time the second voltage has the value Z YX~— up to a con- 
stant. The voltages involved have a full-scale range of 0 to 150 volts and 
accuracy of 1% of full scale is claimed. The “‘settling time’’ of the pulse 
converter is given as 30 milliseconds and presumably this is the limit on the 
repetition rate of the unit. 

F. J. M. 


1124. E. H. Fritze, “Punched card controlled aircraft navigation com- 

puter,” I. R. E., Proc., v. 41, 1953, p. 734-742. 

The computer described is for the purpose of enabling an aircraft pilot 
to fly to an arbitrarily selected point. It provides indications of the heading 
to be flown, the distance to the selected point and the displacement of the 
selected point from a ‘‘master”’ ground station. In addition, the equipment 
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indicates the directions from the aircraft to the master station and an auxili- 
ary station. In flight test, the distance indications obtained during a 200 mile 
flight were in error by less than 23 miles. 

Input data to the computer are obtained from radio receivers which in- 
dicate the directions to a ‘‘master’’ and an auxiliary omnirange station. 
Such stations have been installed throughout the United States by the Civil 
Aeronautics Administration (CAA). The device is also capable of employing 
as input data the range and bearing of a single dme omnirange station. dme 
is a system planned for future installation by the CAA. Through its use, an 
aircraft pilot can determine the distance to a station in addition to the di- 
rection found by means of the omnirange system. 

A substantial amount of auxiliary data is needed by the computer to 
solve its problem. The problem is basically one of triangulation. A unique 
feature is that such auxiliary data is stored on plastic punched cards. By 
positioning the card to reveal the identities of the two stations on which the 
triangulation is based, a set of switches is caused to insert the auxiliary data. 
This includes information which tunes the receivers to the selected stations, 
the rectangular coordinates of one station relative to the other, and a cor- 
rection to account for the difference in magnetic north at the two stations. 

The computer is of conventional analog type employing servomechanisms 
and resolvers. The basic analogy is that of alternating voltage amplitude to 
distance. Sign is indicated by the time phase of the alternating voltage with 
respect to a reference. The most severe limitation on accuracy is the uncer- 
tainty in the bearing data provided by the omnirange system. 

L. H. O'NEILL 
Columbia University 
New York, New York 


1125. R. M. Howe & V. S. HANEMAN, JR., “The solution of partial differen- 
tial equations by difference methods using the electronic differential 
analyzer,” I. R. E., Proc., v. 41, 1953, p. 1497-1508. 

There exist a number of methods for solving linear partial differential 
equations. Frequently one can use separation of variables and reduce the 
problem to ordinary differential equations, which equations can then be 
solved on an electronic differential analyzer (cf. references 1-4 of the paper). 
Another way is to replace the partial differential equation by finite differ- 
ences in all variables and use an arrangement such as a resistance lattice." 
In the present article the authors consider various partial differential equa- 
tions (heat equation in one and two variables, wave equation in one dimen- 
sion, equation of the transverse vibrations of a beam) and replace the spatial 
derivative by finite differences obtaining a differential-difference equation. 
For example, in the case of the heat equation 


C(x) * = 2[xw = | +f(x,t) 
they obtain 
Ca(dun/dt = [Kay3/ (Ax)? ] (agi — Un) —[Ko-4/(Ax)*](un — tn) + fa. 
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A circuit is given solving this equation using operational amplifiers and 
passive elements. Similar differential-difference equations and circuits appear 
for the other equations considered. Sample solutions are given and com- 
parisons with the solutions obtained by the method of separation of variables 
are examined. 
K. S. MILLER 

New York University 

New York, New York 


1F, J. Murray, The Theory of Mathematical Machines. New York, 1947, p. 77. 


1126. R. JENKins, H. W. Brouau, & B. H. Sac, ‘Prediction of tempera- 
ture distribution in turbulent flow. Application of the analog computer,” 
Industrial and Engineering Chemistry, v. 43, 1951, p. 2483-2486. 

A solution obtained on the California Institute of Technology analogue 
computer is compared with a calculated result. The agreement is described as 

2 


1127. E. F. Jonnson, “‘A pneumatic process analog for instruction and re- 
search,”’ Industrial and Engineering Chemistry, v. 43, 1951, p. 2708-2711. 
This instrument is used both for classroom demonstrations and in re- 

search as a check against the mathematical analysis. 


1128. C. A. MENELEY & C. D. Morritt, “Application of electronic differ- 
ential analyzers to engineering problems,” I. R. E., Proc., v. 41, 1953, p. 
1487-1496. 

This is essentially an expository paper. It considers analog electronic 
differential analyzers of the classical type as well as those involving non- 
linear elements. A brief description is given of linear computing elements 
(operational amplifiers used as summers, integrators, etc.) and their basic 
computing circuits ;} nonlinear computing elements such as multipliers and 
function generators; and input-output elements such as recorders and 
plotters. Some typical problems in dynamics are outlined mathematically, 
and closed loop block diagrams for their solution are given. A representative 
list of applications of electronic differential analyzers is included. Finally, a 
few brief remarks are made concerning the amount of equipment required 
for a typical general purpose analyzer. 

K. S. MILLER 

New York University 

New York, N. Y. 


1J. R. RaGazzint, R. H. RANDALL & F. A. RussELt, “Analysis of problems in dynamics 
by electronic circuits,” I. R. E., Proc., v. 35, 1947, p. 444-452. 


1129. WiLLi1AM L. Morris, ‘Analogical computing devices in the petroleum 
industry,” Industrial and Engineering Chemistry, v. 43, 1951, p. 2478- 
2483. 

Three analog devices are described. One is the Phillips ‘‘66 Spectrocom- 
puter,” which is a direct current battery device for solving simultaneous 
linear equations by the Gauss-Seidel method. A second device is a servo de- 
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vice for solving an algebraic equation. The third is an electrical network 
analogue for flow problems. 


1130. D. M. SwIncLe, “Nomograms for the computation of tropospheric 
refractive index,” I. R. E., Proc., v. 41, 1953, p. 385-391. 


1131. P. R. Vance & D. L. Haas, “An input-output unit for analog com- 

puters,” I. R. E., Proc., v. 41, 1953, p. 1483-1486. 

The device described can be used as a recorder in plane cartesian co- 
ordinates or as a curve tracer type of function generator. The unit consists 
of a drum on which is wrapped a piece of graph paper. The drum’s rotation 
corresponds to one variable and the axial movement of a carriage carrying 
either a pen or a potentiometer type transducer corresponds to the other. 
Movements in both coordinates are controlled by servos and static accuracies 
within 0.2 per cent of full scale are claimed. 

For recording purposes the unit functions in the same way as a conven- 
tional plotting board. When used as a function generator the desired func- 
tion is plotted with conducting ink or a soft lead pencil and attached to the 
drum. The transducer, which consists of a printed circuit potentiometer card, 
makes electrical contact with the curve and provides the carriage servo with 
the error voltage required to follow the curve as the drum turns. The output 
is furnished by a potentiometer driven by the carriage motion. Errors in 
following the curve, caused by dynamic limitations of the carriage servo, 
are compensated by adding the servo error voltage to the output potentiom- 
eter voltage. The over-all accuracy seems to be better than one per cent 
within the range for which the unit is intended. The device was designed for 
use with the GEDA computer. 

ROBERT BERNSTEIN 
Columbia University 
New York, New York 


1132. G. B. WALKER, “On the electric field in a multi-grid radio valve,” 
Inst. Elect. Eng., Proc., v. 98, part III, 1951, p. 64-67. 
The use of electrolytic tanks for this purpose is described. 


NOTES 


160.—INVERSE INTERPOLATION FOR THE DERIVATIVE IN THE COMPLEX 
PLANE. In a recent note! formulas were given for finding the argument 
for which a function f(x) has a given derivative f’(x), when that function 
is tabulated for x at equal intervals h. Those formulas are still applic- 
able when dealing with an analytic function f(z) which is tabulated in 
the complex plane, so long as the arguments lie equally spaced upon any 
straight line in the complex plane. But for f(z) tabulated over a Cartesian 
grid x + zy of length h, greater accuracy may be had by locating the argu- 
ments 2 = Zo + kh closer together by choosing k to be small (generally com- 
plex) integers. Thus the problem is to find P, or z = 29 + Ph, when given the 
values of f’, =f’(s) =f’ (20+ Ph) and f,; =f(%) at any conveniently 
located points z. We choose the following configurations of points z for the 
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n-point cases, m = 3(1)7, where the value of m depends upon the number of 
points z; required for direct interpolation: 


3—Point 4—Point 

a a tay 

Zo 21 20 21 

5—Point 6—Point 

2 2145 22: 

Zo 2% 22 i 214i 
20 21 Z2 

7—Point 

225 


Ze Zi4¢  Zeys 
Zo 21 22 


The formula for P in terms of 7, s, ¢, u and 2» is identical with that for p 
given previously,’ namely, 


P=r — rs + r°(2s? — t) + r4(— 55° + Sst — u) 
+ r5(14s4 — 21st + 3 + 6su — 2) 
+ r6(— 4255 + 84s%t — 28st? — 28s’ + 7iu + 7sv) +---. 


But the ¢, s, ¢, « and v are now defined as follows: 
3—Point 
r= {2hf’, + 2(1 —afo— (1 — fi — (1 — Ofi}/(2D), 
s=t=u2zo0=0 
where D = — 2ifo + (1 + af: — (1—d/fi. 
4—Point 
r= {2hf’, + 3(1 — i)fo + 2¢f1 — 2fi — (1 — fii} /(2D) 
s = {3(1 + 4)fo — 3(1 — f1 + 3(1 — Of; — 3(1 + Dfs}/(2D), 


my=ny = 0), 
where D = — Aifo + (3 + if = (3 ™ i)f i + 2th 145. 
5—Point 
r = {10hf’, + 5(4 — 3t)fo + 20if1 — (2 + ife — 4(2 + afi 
— 10f14:}/(10D), 
s = {15(1 + 3i)fo — 15(5 — afi + 6(1 — 20) fe + 3(13 — af; 
+ 15(1 — 34)f14s}/(10D) 
t= {—5(1 + afo+ 10(1 — afi + (1 + S)fe — 2(3 — af; 


+ 10tf14:}/(SD), u=v=0 
where 10D = 5(3 — 11%)fo + 20(3 + 2%)f: — 3(3 — a)fe 


— 4(9 + 24)f; — 10(3 — 24) frsi- 


use 


of 


17! 
no 


sul 


an 
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6—Point 
r = {20hf’, + 40(1 — 4)fo + 16(1 + 2%)f, — (3 — afe — 16(2 + df; 
— 20(1 — afin — (1 — 34) fs} /(20D) 
s = {45(1 + a)fo — 24(2 — 31)f, — 3(2 + Si)f2 + 24(3 — 24)f; 
— 54(1 + a)fi4s — 3(3 + 24)f2x} /(8D) 
t= {— 4fo — 201 + S)fi t+ fe — 201 — HDi t+ Giri + fos} /(D) 
u = {S(1 — t)fo + 4(3 + afi — (1 — 2%)f2 — 4(1 + Sf; — 10(1 — firs 


— (2 — ifes}/(8D) 
v=0 


where 4D = — 30%fo + 32f: — (1 — 3t)f2 — 32f; + 24¢f145 + (1 + fax. 
7—Point 
r = {20hf’, + 4(12 — 114)fo + 40(1 + afi — (1 — 7i)f2 — 8(3 + Ai)f, 
— 20(3 — t)fiss + Sifays — (3 — t)fes}/(20D) 
s = {9(7 + 94)fo — 12(7 — 154)fi — 33f2 + 12(12 — of; 
— 18(3 + 134+) fi4s — 36f244 — 152s} /(8D) 
t= {— (17 + S)fo — 4(4 + %)f. + (5S — Gt)fe — 2(11 — %%)f; 
+ 6(7 + Si)fiss + 6(1 — i) fers + (2 + 3)fes}/(2D) 
u = {(21 — 134)fo + 20(3 + afi + (3 + 144)f2 + 4(3 — 119)f, 
— 1019 — t)fise + 16%f24¢ — 3(2 + @)f2s}/(8D) 
v = {— 3(1 — 3ifo — 12(2 — afi + 3(2 + fe + 6(1 + Sif, 
+ 30(1 — i)fiss — 6(1 + t)fers + 3f2s}/(20D) 
where 20D = 2(8 — 991)fo + 32(9 — 24)f; + (27 + 314)f2 — 96(2 + a)f; 
— 40(4 — 7i) fin, + 32(1 + afar, — (11 — 154) fax. 
H. E. SAaLtzer 
Diamond Ordnance Fuze Labs. 
Washington 25, D. C. 


1H. E. Satzer, “Formulas for finding the argument for which a function has a given 
derivative,” MTAC, v. 5, 1951, p. 213-215. 


161.—TuHE NUMERICAL INTEGRATION OF x(t) = G(x). The formula 
used by EcKERT, BROUWER, & CLEMENCE! in their monumental integration 
of the equations of motion of the five outer planets, is a modification of one 
used by CowELt in his investigation of the motion of Halley’s comet from 
1759 to 1910. In order to get the next position in the ephemeris, use is made 
not only of those already obtained but of first and second summations. 

It seems desirable to have formulas for such work that do not involve 
summations and such are given here. 

Let x = f(t) be a solution of 


x(t) = G(x) 
and put F(t) = G(f(#)). We then have 
f(a + w) = f(a) + wf’ (a) + wf" (@) + ---, 
f(a — w) = f(a) — wf’ (a) + dwf"(a) + ---. 
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We add, express derivatives in terms of forward differences, and then 
replace differences by their values in terms of forward function values, all 
according to Newton’s scheme. Changing w to —w, we see that we can write, 
assuming mth differences of F to be constant, 


f(a + w) = 2f(a) —f(a-—w)+ we, 


where 
Q = AoF(a) + AiF(a — w) + a2F(a — 2w) +---+ AnF(a — nw). 
Let f(t) = #@. Then F(#) = 2, and 
(a + w)? = 2a? — (a — w)*? + w'2(Ao + Ai + Az +---+ An), 


so that 
Ao+Ai+A2+:::+ An = 1. 


Thus if C is the common denominator of the A’s, the sum of their 
numerators is also C, giving a very simple check. 
The following table gives the coefficients for n = 5(1)9. 











n 5 6 7 8 9 
c 240 1 20960 1 20960 36 28800=10! 36 28800=10! 
CAy 317 1 68398 1 76648 55 37111 57 66235 
CA, —266 —1 85844 —2 43594 —92 09188 —112 71304 
CAe 374 3 17946 4 91196 213 90668 296 39132 
CA; —276 —3 11704 —6 00454 —313 23196 —505 69612 
CA, 109 1 84386 4 73136 308 31050 597 00674 
CA, —18 —60852 —2 34102 —203 32636 —492 02260 
CAs 0 8630 66380 86 46188 278 92604 
CA, 0 0 —8250 —21 48868 —103 97332 
CAs 0 0 0 2 37671 22 99787 
CAy 0 0 0 0 —2 29124 
Sum 800 6 79360 12 07360 666 42688 1252 98432 
—560 —5 58400 —10 86400 —630 13888 —1216 69632 
240 1 20960 1 20960 36 28800 36 28800 


The formulas were checked by using simple functions, for instance in the 
case n = 9, the functions f(t) = #° and f(t) = @, with w = 1. 

The formula used by Eckert, Brouwer, and Clemence is for » = 9, and 
by proper transformations it was carried over into a form that avoided sum- 
mation. The coefficients obtained were identical with those directly arrived 
at. Several of the coefficients in the original formula had ten figures in the 
numerators, while none of those above exceeds eight. 

KENNETH P. WILLIAMS 
Indiana University 
Bloomington, Indiana 
1W. J. Eckert, Dirk Brouwer & G. M. CLEMENCE, “Coordinates of the Five Outer 


Planets, 1653-2060,” Astronomical Papers Prepared for the Use of the American Ephemeris 
and Nautical Almanac, v. 12, Washington, 1951. 
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162.—HIsToRICAL NOTE ON Root—FInpinG Macaine. A machine which 
escaped the attention of J. S. FRAME! is one proposed by Frank T. FREE- 
LAND.” The machine is based upon a method for formation of linkages for the 
successive powers, published by Freeland in American Journal of Mathe- 
matics, v. 3, No. 2. 

The machine provides m levers L;, so connected by linkages that a dis- 
placement x of L; produces a displacement x* of Li, k = 2 to nm. On each Ty 
is a pulley P, with axis parallel to and distance $c, from the axis of L;, where 
c, is the coefficient of x*. An inextensible cord is passed around the pulleys 
with the two portions of the cord on P; parallel, perpendicular to L; and on 
the same side of it, using fixed pulleys to change cord direction if necessary. 
To the end of the cord is fixed a pointer brought parallel to a scale marked 
off in units of x. 

The pointer is set at co with x = 0. As L; is displaced through values of x, 
the pointer indicates on the scale the corresponding values of the polynomial. 
A real root is indicated by zero and the real part of a pair of imaginary roots 
by a change of direction of the pointer. 

Another machine is described for quadratic equations based on a linkage 
for x? suggested by SyLvEsTER. A historical sketch is presented in an ap- 
pendix, describing machines proposed by CLArRAuvT in 1820 and DE Roos in 
1879 and mentioning the analytic engine of CHARLES BABBAGE. Finally, 
linkages are described, one for transporting a dimension parallel to itself and 
one for forming a product based on the differences of squares. 

Morris RUBINOFF 
University of Pennsylvania 


Moore School of Electrical Engineering 
Philadelphia, Pa. 


1 J. S. Frame, “Machines for solving algebraic equations,” MTAC, v. 1, 1945, p. 337- 
353. 

? FRANK T. FREELAND, “A machine for the solution of the equation of the m-th Degree,” 
Engineers’ Club of Philadelphia, Proc., v. 2, Feb. 23, 1880. 


QUERIES—REPLIES 


50. INCOMPLETE HANKEL Function. (Q. 43, v. 8, p. 51). 
Put 


z= (x*— 1), s=y(— 1), u = yt = (s?+ y*)! 


to obtain 


fle.) = [7 (@ — 1) 


eo 
= f u—e-“ds 
yz 
° « yz 
= f u— cos uds — if u— sin uds — f u—ds 
0 0 0 


yz ye 
+ f u-(1 — cos u)ds + i f u- sin uds. 
0 0 
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The first two integrals represent Bessel functions, and the third is an ele- 
mentary integral. Thus 


f(x, y) = {4Vo(y) — sinhz + C(y, ys)} — i{4eJo(y) — S(y, ys)}. 


Tables! of the integrals C and S have been reviewed in RMT 651 (MTAC, 
v. 3, 1948-49, p. 479-482). 


A. E. 


1 HARVARD UNIVERSITY, COMPUTATION LABORATORY, Annals, v. 18, 19: Tables of 
Generalized Sine- and Cosine-Integral Functions, Parts I and II, 1949. 


CORRIGENDUM 
V. 4, p. 29, 1. -13, for xx read;11. 








